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Chapter 1 

Introduction and Motivation 



It is widely accepted that the notion of spin was introduced by Uhlenbeck and Goudsmit to explain 
the data on the energy levels of alcaline atoms. The full story is, however, more intriguing and instructive. 
Pauli, in early 1925, decided that the electron has an extra "classically undescribable" quantum number, 
and formulated his famous exclusive principle, according to which two electrons cannot be in the same state. 
This principle made it possible to avoid the difficulties in interpretation of atomic spectra and explain the 
counting of levels in weak and strong magnetic fields. 

Pauli was not, however, willing to make the big jump that the electron has an intrinsic angular momentum 
^h. It was R.L. Kronig, who proposed this idea first. However, this idea was not well received by Pauli, 
as well as in Copenhagen where Kronig went visiting. There was also a problem about the spacing of the 
levels which gave doubts to Kronig himself. Then in the fall of 1925, Uhlenbeck and Goudsmit, in Leiden, 
proposed the same idea which they sent for publication to Naturwissenshaften [||]. After discussions with 
Lorentz, they tried to withdraw their paper, but it was too late (fortunately) and it was published! An 
excellent description of the history of spin and statistics can be found in Ref. Q. 

Presently, spin is a powerful and elegant tool which plays a crucial role in both high energy physics and 
statistical mechanics. Spin is one of the most fundamental properties of elementary particles because it 
determines their symmetry behavior under space-time transformations. 

High energy experiments with polarized beams or final-state spin effects provide often most deep insights 
into the properties of elementary particles and their interactions. For example, the world's best measure- 
ments of the weak mixing angle, sin 6w, have been provided by the SLD experiment at SLAG by using the 
left-right asymmetry in polarized e^e~ scattering, as well as by the LEP results on the forward-backward 
asymmetry for 6-quark final states [^. (For a review, see Refs. ^.) 

It is interesting to note that spin appeared in statistical mechanics also in 1925. The spin model in a 
magnetic field was first solved in one dimension by E. Ising and for that reason it now bears his name. 
In 1944, Onsager first computed the free energy for the two-dimensional Ising model. One of the most 
popular subjects of investigation in the modern statistical physics are critical phenomena in spin systems. 
Present status of our knowledge about the two- and three-dimensional equilibrium spin systems related to 
the Ising, Heisenberg and 0{N) universality classes is discussed in Ref. [^. Recent review on the spin glass 
models can be found in Refs. |, |^]. 

It is well known that there is a close relation between the Quantum Field Theory (QFT) and Statistical 
Mechanics (SM). First, an external similarity there exist: the generating functional of a QFT in the Euclidean 
formulation looks the same as the partition function of corresponding statistical model. This similarity is, 
however, rather formal because both QFT and SM deal with an infinite number of degrees of freedom, and 



further definitions are always needed to remove corresponding divergences. A satisfactory understanding of 
the connection between QFT and SM was reached only when the ideas of the scaling observed in investigation 
of critical behavior of SM models were reconsidered in the general renormalization-group (RG) framework 



by Wilson |11, 12]. Using the field-theoretical methods, it was possible to explain the critical behavior of 
most of the systems and their universal features; for instance, why fluids and uniaxial antiferromagnets 
behave quantitatively in an identical way at the critical point. 

On the other hand, the RG theory of critical phenomena provides the natural framework for defining 
quantum field theories at a nonperturbative level, i.e., beyond perturbation theory (see, e.g., Ref. |jl^). In 
particular, the Euclidean lattice formulation of gauge theories proposed by Wilson |14, 15| provides a nonper- 
turbative definition of Quantum Chromodynamics (QCD), the theory of strong interactions of elementary 
particles. QCD is obtained as the critical zero-temperature (zero-bare-coupling) limit of appropriate four- 
dimensional lattice models and may therefore be considered as a particular four-dimensional universality 
class (see, e.g., Refs. |l^). Wilson's formulation represented a breakthrough in the study of QCD, 

because it lent itself to nonperturbative computations using statistical-mechanics techniques, for instance 
by means of Monte Carlo simulations (see, e.g., Ref. ||l8|| ). 

In this thesis, we study some spin effects in QCD and recurrence lattices with multi-site exchanges. In the 
framework of QCD, we consider the azimuthal asymmetries in heavy flavor production in the lepton-nucleon 
deep inelastic scattering. 

Investigation of the heavy flavor production plays a crucial role in QCD. This is because, for a sufficiently 
heavy quark, the cross sections are calculable as a perturbation series in the running coupling constant a^, 
evaluated at the quark mass. Thus, measurements of the heavy flavor production provide an excellent testing 
ground for perturbative sector of QCD. Moreover, the charm and bottom production is a good probe of the 
structure of the target hadron. In particular, the heavy quark photoproduction is a viable way to measure 
the gluon structure functions (both polarized and unpolarized), while the leptoproduction process is very 
sensitive at large Bjorken x to the intrinsic charm content of the target. 

In the framework of perturbative QCD (pQCD), the basic spin-averaged characteristics of heavy flavor 
hadro-, photo- and electroproduction are known exactly up to the next-to-leading order (NLO). During the 
last fifteen years, these NLO results have been widely used for a phenomenological description of available 
data. At the same time, the key question remains open: How to test the applicability of QCD at fixed 
order to heavy quark production? The problem is twofold. On the one hand, the NLO corrections are large; 
they increase the leading order (LO) predictions for both charm and bottom production cross sections by 
approximately a factor of two. For this reason, one could expect that higher-order corrections, as well as 
nonperturbative contributions, can be essential, especially for the c-quark case. On the other hand, it is 
very difficult to compare pQCD predictions for spin-averaged cross sections with experimental data directly, 
without additional assumptions, because of a high sensitivity of the theoretical calculations to standard 
uncertainties in the input QCD parameters. The total uncertainties associated with the unknown values of 
the heavy quark mass, m, the factorization and renormalization scales, ^ip and /x/?, h.QCD and the parton 
distribution functions are so large that one can only estimate the order of magnitude of the pQCD predictions 



for total cross sections at fixed target energies [19, 2C]. 

At not very high energies, the main reason for large NLO cross sections of heavy flavor production in 
[21, |2^, 7* (7 [^], and gg [24, 25, 26, ^] collisions is the so-called threshold (or soft-gluon) enhancement. 



This strong logarithmic enhancement has universal nature in the perturbation theory since it originates 
from incomplete cancellation of the soft and collinear singularities between the loop and the bremsstrahlung 
contributions. Large leading and next-to-leading threshold logarithms can be resummed to all orders of 



perturbative expansion using the appropriate evolution equations 29, 30 1. Soft gluon resummation of 
the threshold Sudakov logarithms indicates that the higher-order contributions to heavy flavor production 
are also sizeable. (For a review see Refs. |3^, 33|). 

Since production cross sections are not perturbatively stable, it is of special interest to study those 
(spin-dependent) observables that are well-defined in pQCD. A nontrivial example of such an observable 
was proposed in Refs. 35, 36, ^] where the azimuthal cos2ip asymmetry in heavy quark photo- and 
leptoproduction has been analyzed Q. In particular, the Born level results have been considered |3j, ^ and 
the NLO soft-gluon corrections to the basic mechanism, photon-gluon fusion (GF), have been calculated 
| |35| , |36[ |. It was shown that, contrary to the production cross sections, the cos 2^9 asymmetry in heavy 
flavor photo- and leptoproduction is quantitatively well defined in pQCD: the contribution of the dominant 
GF mechanism to the asymmetry is stable, both parametrically and perturbatively. This fact provides the 
motivation for investigation of the photon- (heavy) quark scattering (QS) contribution to the (/^-dependent 
lepton-hadron deep inelastic scattering (DIS). 

In the present thesis, we calculate the azimuthal dependence of the next-to-leading order (NLO) 0(aem"s) 
heavy-flavor-initiated contributions to DIS. To our knowledge, pQCD predictions for the (/^-dependent 'y*Q 
cross sections in the case of arbitrary values of the heavy quark mass m and are not available in the 
literature. Moreover, there is a confusion among the existing results for azimuth-independent 7*Q cross 
sections. 

The NLO corrections to the (^-independent lepton-quark DIS have been calculated (for the first time) a 
long time ago in Ref . |3^] , and have been re-calculated recently in |3^] . The authors of Ref. ||3^ conclude 
that there are errors in the NLO expression for u'-^-* given in Ref. [38| 0. We disagree with this conclusion. 
It will be shown below that a correct interpretation of the notations for the production cross sections used 



in |38| leads to a complete agreement between the results presented in Refs. and present thesis. 

As to the (/^-dependent j*Q cross sections, our main result can be formulated as follows. Contrary to the 
basic GF component, the QS mechanism is practically cos 2(/?- independent. This is due to the fact that the 
QS contribution to the cos 2ip asymmetry is absent (for the kinematic reason) at LO and is negligibly small 
(of the order of 1%) at NLO. This fact indicates that the azimuthal distributions in charm leptoproduction 
could be a good probe of the charm density in the proton. 

Then we investigate the possibility of measuring the nonperturbative intrinsic charm (IC) |^ using the 
cos 2(p asymmetry. Our NLO analysis of the hadron level predictions shows that the contributions of both 
GF and IC components to the cos 2(p asymmetry in charm leptoproduction are quantitatively well defined: 
they are stable, both parametrically and perturbatively, and insensitive (at > rm?) to the gluon transverse 
motion in the proton. At large Bjorken x, the cos 2ip asymmetry could be a sensitive probe of the intrinsic 
charm content of the proton. 

We have also considered the contribution to azimuthal distributions of the perturbative charm density 
within the variable fiavor number scheme (VFNS) |Q, 44] ^. Main result of our analysis is that the charm 



densities of the recent CTEQ [^5|] and MRST [^6[ sets of parton distributions have a dramatic impact on 
the cos 2(p asymmetry in the whole region of x and, for this reason, can easily be measured. 



^The well-known examples are the shapes of differential cross sections of heavy flavor production which are sufficiently stable 
under radiative corrections. 

^For more details see PhD thesis ji^], pp. 158-160. 
The notion of the IC content of the proton has been introduced over 25 years ago in Refs. 0,0]. This nonperturbative 
five-quark component, \uudcc), can be generated by gg cc fiuctuations inside the proton. 

*The VFNS is an approach alternative to the traditional fixed flavor number scheme (FFNS) where only light degrees of 
freedom (u, d, s and g) are considered as active. Within the VFNS, potentially large mass logarithms are resummed through 
the all orders into a heavy quark density which evolves with according to the standard evolution equation. 



Concerning the experimental aspects, azimuthal asymmetries in charm leptoproduction can, in principle, 
be measured in the COMPASS experiment at CERN, as well as in future studies at the proposed eRHIC 
1 4^, Hi and LHeC H] colhders at BNL and CERN, correspondingly. 

Another topic of our thesis are critical phenomena in spin systems defined on the recurrence lattices 
with multi-site exchanges. It is well established that the thermodynamic properties of a physical system 
can be derived from a knowledge of the partition function. Since the discovery of statistical mechanics, it 
has been a central theme to understand the mechanism how the analytic partition function for a finite-size 
system acquires a singularity in the thermodynamic limit when the system undergoes a phase transition. 
The answer to this question was given in 1952 by Lee and Yang in their seminal papers ||5^, It was 
shown that phase transitions occur in the equilibrium systems in which the continuous distribution of zeros 
of the partition function intersects the real axis in the thermodynamic limit. For anti- ferromagnetic Potts 
models, by contrast, there are some tantalizing conjectures concerning the critical loci, but many aspects 



still remain obscure |53, 53, pM. Recently, the Yang-Lee formalism has also been applied for investigation 



of nonequilibrium phase transitions [^ ]. 

Presently, the investigation of the partition function zeros is a powerful tool for studying phase transition 
and critical phenomena. Particularly, much attention is being attached to the study of zeros of partition 
function of helix-coil transition of biological macromolecules [56, 57, 58, p9|] . 

In this thesis, we investigate helix-coil phase transition for polypeptides and proteins in thermodynamic 
limit on recursive zigzag ladder with three-spin interaction. We use recursive lattices because for the models 
formulated on them the exact recurrence relations for branches of the partition function can be derived. For 
classical hydrogen bond {N — H ■ ■ ■ O = C), we have got the Yang-Lee zeros corresponding to helix-coil phase 
transitions for polypeptides and proteins in thermodynamic limit. We also take into account a non-classical 
helix-stabilizing term describing a hydrogen bond of the type Ca — H ■ ■ ■ O. For this case we obtain folding 
and quasi unfolding of the order parameter (degree of helicity 0). Applying multi-dimensional mapping 
on zigzag ladder, we got Arnold tongues for non-classical helix-coil phase transition for neutral points of 
mapping 0. 

There are two types of modulated phases: commensurate and incommensurate ones. For commensurate 
phases, when (f = ^tt, the so-called Arnold tongues there exist. Typically, for multi-dimensional maps, the 



border of such regions (Arnold tongues |6^) splits into two branches in the parameter space. 

Our main result is that we get two qualitatively different behaviors for the degree of helicity that 
depend on input parameters. The first regime presents a low-temperature helix structure which melts at 
higher temperatures. We observe that the presence of a non-classical (Ki) interaction sensibly enhances the 
melting temperature, and the transition is smooth. In second case, the presence of non-classical interaction 
leads to a remarkably different low-temperature behavior. In this regime, an quasi unfolding transition takes 
place for T ^ as well, like to cold denaturation [31|. We point out that our results are meaningful for 
long chains since, for such chains, a thermodynamic limit is involved. Note that unfolding of biopolymer 
has also been observed in phenomenological model [p2|, Monte Carlo simulation l63], Bethe approximation 



1 64], and for a short chain in Distance Constraint Model [65|. 

In this thesis, we also investigate magnetic properties of the ^He. Investigation of magnetic phenomena 
and magnetic properties of materials has a long history |^6[ . The theory of magnetism and related problems 
composes a fast and rather advanced field of research in the modern theory of condensed matter physics 
intimately linked to many other fields of physics, mathematics, biophysics, chemistry and materials science. 
These investigations have a wide application in various fields of electronics, computer techniques e.t.c. The 



^Neutral points are defined by condition that eigenvalues of tfie mapping Jacobian, A, lie on the unit circle, A = e^' 



unexpected discovery of cooper-oxide high-Tc superconductors in 1986 |67] not only aroused hopes that one 
day we wih have at our disposal materials which exhibit superconductivity at room temperature, but also 
opened a new stage in the studies of magnetic phenomena. This is because there is a strict evidence that 
two-dimensional anti-ferromagnetism is one of the key components of high-temperature superconductivity. 
One of the most remarkable achievements in this field is the progress in the studies of magnetism in soHd 
^He. 

Solid and fluid '^He films absorbed on the surface of graphite have attracted extensive attention since 
(at low temperatures and high pressures, due to the light mass of helium atoms) it is a typical example 
of a two-dimensional frustrated quantum-spin system p8| , |6S| ]. The nuclei of ^He are fermions with spin 
1/2. It's reasonable to assume that solid '^He is a system of localized identical fermions. The microscopic 
theory of magnetism for such systems is based on the concept of the permutation of particles. In the 
films under consideration, the nuclei of '^He form a system of quantum 1/2 spins on a triangular lattice. 
We know experimentally that the three-particle interactions dominate in this system. Transition from 
ferromagnetic behavior to antiferromagnetic one takes place when the coverage (density) of ^He atoms 
decreases. The explanation is suggested in terms of multiple— spin exchanges (MSE). In a dense clode- 
packed solid, ferromagnetic three-spin exchange is dominant |7[1[ |. At lower densities, ferromagnetic three- 
and-five spin exchanges compete with antiferromagnetic four-and-six spin exchanges and lead to a frustrated 
antiferromagnetic system. The MSE produce frustration by themselves and a strong competition between 
odd— and even— particle exchanges is also responsible for the frustration For description of solid ^He 
films, one can use the dynamical system approach with MSE model that leads to appearance of various 
ordered phases and magnetization plateaus and one period doubling |72, |7^. The study of the above 
mentioned magnetization plateau is one of the main directions of present-day activity in the field of non- 
trivial quantum effects in condensed matter physics. Despite the purely quantum origin of this phenomenon, 
it was shown recently that magnetization plateau can appear in the Ising spin systems as well exhibiting in 
some cases fully qualitative correspondence with its Heisenberg counterpart ||7^, |7|, 77 1. 

Using the dynamical system approach with MSE on the recurrent lattices (square, Husimi, hexagon), 
we obtain magnetization curves with plateaus (at m = 0,m = 1/2, m = 1/3 and m = 2/3) and one period 
doubling. 

The next issue of our investigation is the so-called face-cubic model. We have considered a spin model 
with cubic symmetry defined on the Bethe lattice and containing both linear and quadratic spin-spin in- 
teractions. An expression for the free energy per spin in the thermodynamic limit was obtained. We have 
applied the methods of the dynamical systems theory or, more precisely, the theory of discrete mappings. 
In this technique, one exploits the self-similarity of the Bethe lattice and establishes a connection between 
the thermodynamic quantities defined for the lattices with different number of sites. 

We have identified the different thermodynamic phases of the system (disordered, partially ordered and 
completely ordered) in the ferromagnetic case {J > 0, K > 0) with different types of the fixed points of 
recurrent relation. Then we have obtained the phase diagrams of the model which are found to be different 
for Q < 2 and Q > 2. The case of Q < 2 contains three tricritical points while only one tricritical and one 
triple points there exist at Q > 2. 

Our results on the critical phenomena in spin systems defined on recurrence lattices with multi-site 



exchanges are published in Refs. [[Tg, |8lj, |8^. Our studies of the spin effects in QCD are presented 

in Refs. |8|, ||. 

The thesis is organized as follows. In Chapter the multi-dimensional mapping is used for non-classical 
helix-coil phase transition of anti-ferromagnetic Potts model for biopolymer formulated on the recursive 



zigzag ladder. Two qualitatively different behaviors for the degree of helicity are obtained. 

In Chapter ^ three types of the recurrent lattices with MSE are considered as approximation to solid 
^He films. Using methods of the dynamical systems theory, we've got magnetization plateaus, bifurcation 
points, one period doubling behavior and modulated phases at sufficiently high temperatures. 

In Chapter ^, we derive the system of recurrent relations for the face-cubic model on the Bethe lattice. 
We identify different types of the fixed points of the system of recurrent relations with different physical 
phases. 

In Chapter ^, we analyze the QS and GF parton level predictions for the (/^-dependent charm leptopro- 
duction in the single-particle inclusive kinematics. Hadron level predictions for azimuthal asymmetry are 
obtained. We consider the IC contributions to the asymmetry within the FFNS and VFNS in a wide region 
of X and Q^. 

Main observations and conclusions of this thesis are discussed in Conclusion. 



Chapter 2 

Advantage of Recursive Lattices i 



The advantage of recursive lattices is that for the models formulated on them the exact recurrence 
relations for branches of the partition function can be derived. Let us consider the recursive lattices which 
are connected through the sites. As the first example of recursive lattice is an usual chain. One can receive 
the exact recursion relation for the partition function for Ising model. We divide a chain on two equal parts. 
The partition function can be written as follows: 

Z = ^exp(/?/icTo)-(7'(^^o) (2.1) 

where ao is the center of the chain and gn{<^o) is the contribution of each chain branch. gn{(^o) can be 
expressed trough g„_i(c7i), that is, the contribution of the same branch containing n-1 generations starting 
from the site belonging to the first generation: 

5n(o"o) = ^exp(JcJocri + hai)[gn-i{cri)], (2.2) 

where cTj takes values it 1, J is interaction constant and h is the external magnetic field. We introduce the 
following variable 

x„ = ^, (2.3) 

9n (-) 

where we denote gn{<^o) by gn (+) if the spin do takes the value +1 and by gn (— ) if the spin takes the 
value -1. For x„ we can then obtain the recursion relation: 

x„ = /(x„_i). (2.4) 

/(x) is a ratio of two polynomials. We obtain one dimensional dynamic rational mapping. We can get the 
magnetization of a central site through Xn- 

Another example of recursive lattice is the Bethe one. Let us regard the Potts model on this lattice with 
7 coordination number. The Hamiltonian can be written as: 

n = -J 5{<Ji,aj)-HY,K^^^)■ (2.5) 

<i,j> i 

where 5{ai,aj) = 1 for ai = and otherwise, ctj takes the values 1,2,...Q, the first sum is over the nearest- 
neighbor sites, and the second sum is simply over all sites on the lattice. We use the notation K=J/kT 
and h=H/kT. Cutting apart the Bette recursive lattice at the central point we get 7 identical branches. As 



The results considered in this chapter are pubhshed in Refs. |7|, ||. 



usual we can receive one dimensional dynamic rational mapping for partition function. The same ideas can 
be used as for the recursive chain. Denoting gn{co) the contribution of each lattice branch one can receive 
the recursive dynamic relation. Introducing the notation 



.„ = ?^(^ (2,6) 

9n (C^ = 1) 



one can obtain the Potts-Bethe map 



The magnetization of the central site for the Bethe lattice can be written as 

e^^ + {Q~- l)x} 

The situation changes drastically for Q < 2 with antiferromagnetic interactions. The plot of the M (magne- 
tization) versus h (external magnetic field) has a bifurcation point and chaotic behavior at low temperatures 



Mn =< S{ao, 1) >= , ^ (21 



An other example is Husimi lattice. It can be regarded as recursive lattice. The three-site antiferromag- 
netic Ising model on Husimi lattice is investigated in an external magnetic field using the dynamic system 
approach. Making the same procedure for Husimi recursive lattice one can obtain one dimensional rational 
recursive relation for partition function. The full bifurcation diagram, including chaos, of the magnetization 
was exhibited. It is shown that this system displays in the chaotic region a phase transition at a positive 
"temperatures" whereas in a class of maps close to x 4x(l — x), the phase transitions occure at negative 
"temperatures". The Frobenius-Peron recursion equation was numerically solved and the density of the 
invariant measure was obtained l87|l. 



The ladders [g^, |8^, Q also can be regarded as a recursive lattice. They are connected through the bonds 
and have multi dimensional rational mapping for partition function. A zigzag ladder with axial next-nearest- 
neighbor Ising model has attracted many investigators on account of the fact that it is a particularly simple 
model exhibiting quasi specially modulated phases that can be either commensurate or incommensurate 
with the underlying lattice Using the dynamic approach one can receive three dimensional rational 

mapping for partition function. 

2.1 Arnold Tongues in Ising and Potts Models 

Let us regard the anti-ferromagnetic Ising and Potts models on the recursive Bethe lattice connected 
through sites. For Ising model the partition function can be written as: 

Z= ^exp{/iaoK(ao), (2.9) 
where ctq is the central spin, gn{<yo) - the contribution of each lattice branch, h - magnetic field and q - 



coordination number|l92|. gn{c^o) is obviously expressed through gn^i{ai): 

gM = E'M ^^' (2.10) 

0-1 

for q = 3 and interaction between the spins is constant J = — 1. Introducing the notation 

Xn = ^ (2.11) 

9n{-) 




the recursion relation (plo| ) can be rewritten in the form 

Xn = f{Xn~l,T, h). 



(2.12) 



As is known, if the derivative of f{x,T,h) is equal to —1 we have a bifurcation point, corresponding to 

_ 2_ 

the second order phase transition for anti-ferromagnetic model. We defined v = e t and after a simple 
calculation we have got the following system of equations: 



yh + l ^2 



(2.13) 



2vx — 2x^ 



-1 



Eliminating x we obtain the following equation: 

Solving this equation we get: 

- f = -3 In 2 + I + In 1 1 - 6^2 - 3?;^ ± v/(l - 6v^ - Sv^)"^ - 6Av<^^ 



(2.14) 



(2.15) 



This equation define Arnold tongues between paramagnetic and modulated phases with winding number 
w = 1/2. 

The Arnold tongue begins at the temperature of T = when the external magnetic field h = 0, and 
ends (T = 0) at /i = ±3 (see figure 2.1). 

The same procedure we can perform for anti-ferromagnetic Potts model on recursive Bethe lattice with 
Hamiltonian 

-f5H = - S{ai,aj) + hY,H^iA), (2-16) 

<i,j> i 

where ai takes the values 1,2,3. Introducing the notation 

_ 9n(.*) 



5n(l)' 



(2.17) 




0.1 0.2 0.3 0.4 0.5 0.6 

T 

Figure 2.2: Arnold tongue for anti-ferromagnetic Potts model on recursive Bethe lattice with coordination number 
q=3. 



where g'n(l) is the branch of partition function with central spin a = 1 and is the branch of partition 

function with central spin a ^ 1. For the coordination number of the Bethe lattice g = 3 we obtain following 
system of equations 



+ {z + l)x' 
-h+l _^ 



(2.18) 



2(z + l)x - 4x^ 



-1 



here again the derivative of /(x, T, h) is equal to —1 which corresponds to the second order phase transition 
of anti-ferromagnetic model and where z = e~T . The Arnold tongue begins at z = g(\/T7 — 3), when 



external magnetic field h = 1.5, and ends (T = 0) at /i = and h = 3 (see figure 2.2). 



2.2 Multi-dimensional Mapping for the Biological Macromolecules 



The structure of a protein is completely encoded in the amino-acid sequence [p3| . Understanding of the 
folding and unfolding processes of proteins (hetero-polymers) and polypeptides (homo-polymers) is one of 
the current challenges in molecular biophysics. A lot of effort has been devoted to clear up the mystery of 



protein or polypeptide folding nature by using lattice models [94|. These simple lattice models single out 
the formation of a helix structure in protein as the basic mechanism to be understood. Thermodynamics 
of homo and hetero-polymers folding has been investigated in this perspective by introducing a variety of 



different lattices: chains [95, 96, 97|, square lattices [98, 99[, and cubic ones [IOC, 101[. Off-lattice models 

Chaotic behavior in off-lattice 



have been discussed by Irback et al. and Klimov and Thirumalai [102, |10S 
models of hetero-polymers (proteins) and folding and unfolding have been analyzed in two-dimensional 



systems by means of Monte Carlo simulations [ 104 [. The theory of finite-size scaling of helix-coil transition 
was studied by Okamoto and Hansmann [|5^, 57, 58, |59|[ by multi-canonical simulation. They have chosen 



three types of polypeptides with aliphatic neutral amino acids (alanine, valine, and glycine). It was shown 
that a-helix formation in short peptide systems agrees with experimental results |105([ . But proteins are 
composed of different types of monomers. Hydrophobic monomers, such as leucine or proline, try to hide 
their surfaces from the solvent. The simplest protein theoretical model divides the amino acids into two 
categories: hydrophobic (H) and polar (P) surrounded by the solvent [IOC, 107[. Kamtekar et al. [108[ made 



experiments with a variety of hydrophobic (H) and polar (P) amino acids in hetero-polymers and showed 
that a simple code of polar (P) and nonpolar hydrophobic (H) residues arranged in an appropriate order 
could drive polypeptide chains to collapse into globular a-helical folds. By using a simple HP lattice model 
T0| a theory explained the experimental phenomenon of cold denaturation (unfolding) on real proteins 



109 



61]. The study of relaxation processes in biopolymer is of particular significance since the functional abilities 



of thes molecules are related to the dynamical properties 111|. 

We point out that different theoretical models were proposed to study both unfolding and folding of 
proteins p^ . 

From a statistical mechanics perspective different approaches can be attempted to investigate the nature 
of the helix-coil phase transition: from the analysis of Yang-Lee zeroes 
Carlo simulation for finite samples |5q, p% 



51 1, to multicanonical Monte 



II, |59|, ^, |99|, |100|, |10J, |109|, |11[ 

In this thesis we study in thermodynamic limit a model for the helix structure of proteins and polypep- 
tides, where we take into account both the classical hydrogen bond |112] between three a-carbons by using 
CO and NH H-bond connection, and the non-classical H-bonds [|113|| in every Ca — H. The classical (a— 
helix) H-bond is formed in the following way: three neighboring angle pairs [Ca{^i,'4^i),Ca{^i+i,ipi+i) and 
Ca{fi+2,ipi+2)] form a H-bond when rotations are such that the distance between H [N{i — 1) — H] and O 
[O = C{i + 3)] becomes less than 2 A (fig.l). The hydrogen bond is a unique phenomenon in structural 
chemistry and biology. Its functional importance stems from both thermodynamic and kinetic reasons. In 
super molecular chemistry, the hydrogen bound is able to control and direct the structures of molecular 
assemblies because it is sufficiently strong and sufficiently directional. The subject of hydrogen bonding is 
of major interest and remains relevant with each new phase in the kaleidoscope of chemical and biological 
research (see references in |114| ]). Non-classical C — H . . .0 bonds have been recognized to play an important 



role in biological macromolecules (see [115]), and they were for instance observed between water and amino 
acids alanine [Ca — H ■ ■ ■ Oi?2] [ |116[ or between two helices of collagen [117], 

Traditionally, the transition from random coiled conformation to the helical state in DNA, RNA or 
proteins are described in the framework of Zimm-Bragg ]|118f] type Ising model. But this type of one- 



dimensional model cannot account for non-trivial topology of hydrogen bonds ]112]. 



N 

I 

H 



R H 
\ / 



c 

II 





H 

I 

N 



H-bond 



R H 

\ / 



H 

I 

N 



R H 



N 

I 

H 



C 

II 







II 

c 




R H 



N 

I 

H 



Figure 2.3: The backbone of polypeptide or protein. The classical H-bond interaction between N-H and C=0 is 
pointed out by dashed line. 



We show the backbone chain of the polypeptide molecule in fig. 2.3. R (amino acid residue) denotes the 
side chain. Because of the planar structure of the amide group, almost the whole conformational flexibility 
of the polypeptide backbone chain is determined by the rotation angles around the single bonds — Ca 
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Figure 2.4: The zigzag ladder. 3-site Potts H-bond interaction is marked by solid line. 



and Cq, — C which are usually denoted as 99 and i/' respectively. 

We now formulate our model: first of all the angle pairs {ipijtpi) are discretized |119|, the possible Q 
values are labeled by a discrete variable Sj.When three successive rotation pairs (spins) are zero (sj = Sj+i = 
= 0), an H-bond appears which leads to some energy gain. When one of three neighboring spins is not 
zero(sj = *) an interaction with solvent is taken into account. This leads to a three-site interaction Potts 
model [|12C1|] on a zigzag ladder (see fig. |2.4|) The Hamiltonian of the system is written as 



H = -jY,S{si-i,0)S{si,0)6{si+i,0) (2.19) 
A' 



i 

where J is the energy of hydrogen bond, K is the energy of protein-solvent hydrogen bond, Si denotes the 
Potts variable at the site i and takes the values 0, 1, 2, • • • ,Q — 1, Ki is the energy of non-classical H-bond, 



and Aj label each triangle in Fig. 2.4. 

The model we thus introduced is indeed quite a simplified one, but it allows to discuss how non-classical 
bonds compete with classical hydrogen interaction in an idealized setting. We will take advantage of the 
recursive nature of the zigzag ladder: this makes it possible to derive exact recursion relations for branches of 
the partition functions, and in this way statistical properties in the thermodynamic limit may be investigated 
by dynamical systems techniques 1 75 ] . In their simplest realization recursive relations yield one dimensional 
mappings |73| . 

By cutting the zigzag ladder in the central triangle (s_i, sq; ■^i) one gets the partition function associated 
to the hamiltonian ( 2. IS] ) 

E [e-^Z(^\s.i,so)Z(^\so,si)], (2.20) 

{s-i,so,si} 

where 

Ho = -{J - K)5{s^i,0)5{so, 0)5{su0) - Ki (5(s_i, 0) + 5{so, 0) + <5(si, 0)) , (2.21) 

T is temperature (room temperature is T = 0.6^^^), s_i,so,si are spins of central triangle, Z(")(s_i,so) 
and Z(")(so, si) are the parts of partition function corresponding to two branches, n is generation of recursive 
lattice (see Fig. p.4D . By introducing the following notation 

(0,0)= ; (0, *) = ) ; 



Z('^)(*,0) = Zf); Z('^)(*,*) = ZfK (2.22) 



and 



J-K Ki , ^ 

7 = exp ^ ; 2; = exp— , (2.23) 



( ^.20 ) can be rewritten as: 

Z ~ 7^'[^|"¥ + 2(Q - l)z'zt^Z^^^ + {Q- lfz[Z^^^f + {Q- \)z\zff (2.24) 

+2(Q - l)2zZ(")zf ) + (Q - \f{zf\^ 

By applying the "cutting" procedure to an nth generation branch one can derive the recurrence relations 

r r7(n) ry{ri) ry{n) r7(n) 

lor Zj^ , , z-g , Z/4 , 

Zi"^ = 7zZ;"-') + ) = zZf + (2.25) 

3 — ZZ/j^ -\- , Z^ — ^Z/g -|- Z/^ 

If we notice that Zg"^ = z|"^ , and introduce the notation 

^(n) ^(n) 

Z4 Z4 



we can obtain a two-dimensional mapping from ( |2.25| ), 

Xn = fl{Xn~l,yn~l), h{x,y) = ^^^^^ (2.27) 

zy + Q - I 

yn = f2{Xn-l,yn-l), f2{x,y) = "|" ^ 1. 

-zy + Q - 1 

2.3 Helicity and Arnold Tongues for the Macromolecules 

The (dimensionless) order parameter or helicity defined as 

where = (5(s_i, 0)(5(so5 0)(5(si, 0)) (when recursion relations admit a stable fixed point the order parameter 
is independent of the triangle we consider). Since our procedure implicitly involves a thermodynamical limit, 
its biological significance is motivated by the existence of long chains of proteins, like collagen, that may 
exist in the form of three intertwined peptide chains, each containing a thousand of amino acids (we also 
remark that the importance of non-classical H bonds in collagen has been pointed out in [|117[| ). 

Thus our task is that of investigating the asymptotic behavior of recursion relations ( 2.27| ), this allow 



to characterize the macroscopic order parameter as a function of the physical parameters T, J (energy of 
classical H bond), K (energy of protein-solvent H bond) and Ki (energy of non- classical H bond). 

We observe that in the whole range of the parameters triplet {J,K,Ki) the recursion relations ( |2.27| ) 
admit a single (real) fixed point {x,y). An investigation of the Jacobian matrix of the transformation at 
(x, y) moreover indicates that such a fixed point is always stable: thus we do not get any phase diagram 
marked by the stability border for the fixed point like in mean field Ising models with competing interactions 



on hierarchical lattices 121 |. We also point out that we did not observe in our tests any other dynamically 



relevant attracting structure: under iteration of ( 2.27 ) generic initial conditions collapse to the stable fixed 



point. (We mention that for physical values of microscopic parameters there exist complex fixed points 
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Figure 2.5: The order parameter as a function of temperature (T) for different values (Q, J, K,Ki). 



(x, y) at which the absolute value of eigenvalue of Jacobian equals to one. In this case the order parameter 
(0) would be complex too). 



Once {x,y) is determined, the degree of helicity ( 2.28| ) can be computed. Our main result is that we get 
two possible qualitatively different behaviors. The first regime presents a low-temperature helix structure, 
which melts at higher temperatures (see fig j2.5| ), a qualitative feature that may be observed if the non 



classical interaction is absent. By looking at two of the curves in fig j2.5| we observe that the presence of 
a non-classical (i^i) interaction sensibly enhances the melting temperature, and that, coherently with the 
dynamical analysis of recursion relation, the transition is smooth. For other parameter values, the presence 
of non-classical interaction leads to a remarkably different low temperature behavior, with an quasi unfolding 
transition also for T — > 0, akin to cold denaturation see fig.^!^. Real unfolding behavior is when order 
parameter's peak is near one. In fig.^]^ for parameters: Q = 4, J = 0.8, K = 2 and Ki = 3.5 we have 
a strange situation. At low temperature we have a coil. The protein changes conformation upon heating. 
Some of [N{i — 1) — H] and = C(i + 3) in average are near and they form H-bond. In these parameters 
the degree of helicity (order parameter) Q becomes larger. We call that quasi helix. At higher temperature 
the protein becomes coil again. 
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Figure 2.6: Starting from low temperature and upon heating, the protein changes conformation from coil to quasi 
helix, and then at still higher temperature becomes coil (disordered). 



Using the theory of dynamical systems for two-dimensional mapping, we have obtained the separating 
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Figure 2.7: The line separating coil (paramagnetic, disordered) and helix (modulated, ordered)phases. 




Figure 2.8: Arnold tongue with winding number w — j2tV — f"" ^^id Q = 50 for non-classical helix-stabilizing 
interaction. 

line, which divides the coil (paramagnetic, disordered) phase from helix (modulated, ordered) one (see figure 
|2.7|).Two example of Arnold tongues for non-classical helix-stabilizing interaction with Q = 50 for = Stt 



w = To ^"-^ ^ ~ 1^ = i Eire shown on figures |2.8| , p.S| . We point out that our result are meaningful 




Figure 2.9: Arnold tongue with winding number w — |,(^ = |7r and Q = 50 for non-classical helix-stabilizing 
interaction. 



for long chains, since a thermodynamic limit in the statistical model is involved. We notice however that 
unfolding of biopolymer has been observed in phenomenological model |^^, Monte Carlo simulation [|6^ ], 
Bethe approximation |6^, and for a short chain in Distance Constraint Model p5| . 



2.4 Yang-Lee Zeroes for the Biological Macromolecules 



When we take into account only classical hydrogen bound the Hamiltonian of the system is written as 
-/3H = jY^ <5(s,_i, 0),5(si, 0)5(s,+i, 0) + ^[1 - 5{si.i,0)6isi, 0)6{si+i,0)], (2.29) 



We obtain again the two dimensional rational mapping relation for Xn and yn 

Xn = fl{Xn-l,yn-l) 
Vn = f2iXn-l,yn-l), 



(2.30) 



where 



Mx,y) -- 
f2{x,y) 



lx + {Q- 1) 

y + {Q-i) 

x+{Q-l) 

y + {Q-i)' 



(2.31) 



In case of multi-dimensional rational mapping the fixed point x*,y* is attracting when the eigenvalues of 
Jacobian |A| < 1, repelling, when |A| > 1, and neutral, when |A| = 1. So the system undergoes a phase 
transition when 



W - exp {t<f) 
dx 



dy 

^ - exp (up) 







(2.32) 



After eliminating x and y from ( |2.31] ) and ( 2.32| ) we obtain the following equation for the partition 
function zeros 



where 



bQ + bi cos{ip) + 62 cos^((/3) + 63 cos^((/3)[cos(3(/9) + ism{3ip)] = 0, 



^'o = Q'(7-l) + (7-l)7' + Q'(7 + l)+Q7(7'+7-2) 

61 = 2{Q'^ + 2Q{j-l)-Q\j^-l)+^{j^ + j-2)} 

62 = -4(Q-l)(7-l)(2 + Q + 7) 

63 = -8(Q-1)(7-1). 



(2.33) 



(2.34) 



One can solve ( p. 33 ) for fi and find the Yang- Lee zeros of partition function with different parameters 
Q, J and T like in Refs.[|T2^. These parameters are different for each polypeptides and proteins. After 
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Figure 2.10: The Yang-Lee zeros: a) Q=9, J=2.0, T=0.56, K=0.76; b) Q=30, J=2.1, T=0.54, K=0.26; c) Q=40, 
J=2.4, T=0.56, K=0.334; d) Q=50, J=2.5, T=0.6, K=0.1523. 



making discrete values of Q and comparing with Ramachardan and Shceraga[ll£] we confirm that the 
circle in classical helix-coil transition does not cut the real axis. So we have not a real phase transition in 
polypeptides (proteins). According to phenomenological theory of Zimm-Bragg or Lifson-Roig there is only 
pseudo phase transition in 2-site (Ising) model. Our results describe the microscopic theory of helix-coil 
transition of polypeptides or proteins with non-trivial topology of hydrogen bonds and find Yang-Lee zeros of 
pseudo phase transition. Yang-Lee zeroes of helix-coil transition for polyalanine, polyvaline and polyglysine 
was regarded [^]. The authors made Monte Carlo simulation technique and considered polypeptide chain 
up to = 30 monomers and determine the (pseudo) critical temperatures of the helix-coli transition in 
all-atom model of polypeptides. 



Chapter 3 

Fluid and Solid ^He i 



As mentioned in introduction, fluid and solid ^He films absorbed on the surface of graphite have attracted 
extensive attention, since it is a typical example of a two-dimensional frustrated quantum-spin system 
p^ , 39 1 . The first and second layer of the nuclei of ^He forms a system of quantum one-half spins on a 



triangular lattice. The third layer forms a Kagome one |123f| . Many experimental [124| and theoretical 
1 125 1 studies suggest that the exchange of more than two particles are dominated in these systems. For 
such systems a change from ferromagnetic behavior to anti-ferromagnetic takes place. Spin ladder anti- 
ferromagnets have been attracting extensive interest because they have a spin gap. A special type of 
frustration due to cyclic exchange interactions was recently found to be important in the spin ladder material 
LaxCai4_a;Cu2404i |8^. It is experimentally also known that a many-body exchange interaction cannot be 
neglected especially in ^He on graphite [126]. 

Last decade the investigation of magnetization plateaus in a strong magnetic field has taken on special 
significance. The magnetization plateaus are famous for the fact that they are an example of essentially 
macroscopic quantum phenomenon. For the first time, Hida has theoretically predicted an appearance of the 
magnetization plateau for the ferromagnetic-ferromagnetic-antiferromagnetic Heisenberg chain of 3CuCl2-2 
dioxane compound, which consist of the antiferromagnetic coupled ferromagnetic trimers |127]. The values of 
magnetization at which the plateaus appear are quantized to fraction values of the saturation magnetization. 



The theoretical explanation of this fact was been given in 1997 by Oshikawa, Yamanaka and Affleck |12S]. 
These magnetization plateaus were observed as a simple origin in the Ising limit |129]. Geometric frustrated 
quantum magnets are a class of magnetic materials with various unusual properties at low temperature and 
high pressure. Due to strong frustration and quantum effects, these materials may be in principle considered 
as a source of new strongly correlated physics. The most of studied geometric frustrated quantum magnets 
are the Kagome and pyrochlore lattices of antiferromagnetic coupled nearest neighbor spins. As mentioned 
above, the third layer of the nuclei of ^He forms a Kagome lattice. Usually, the antiferromagnetic Kagome 
nets are investigated using numerical simulations [130|. We propose a dynamic approach based on exact 
recursive relations for partition functions. Our method makes possible to research magnetization plateaus, 
bifurcation points and period doubling in anti-ferromagnetic case at low temperatures and high pressures. 



^ The results considered in this chapter are pubhshed in Refs. 



3.1 Ising Model Approach to the Solid ^He System on the Square Re- 
cursive Lattice 



The most general expression for the Hamiltonian with multi spin-exchanges on a triangular lattice is 

H = Hph + 'Hex + T~^z- (3.1) 

The term Hph describes the phonon contribution and is not essential. Hex responses for two-, three-, 
and four exchange interactions. Hz term is responsible for magnetism in solid ^He and is given by Zeeman 
Hamiltonian 

Hz = -J2lhH-ai (3.2) 

i 

where 7 is gyromagnetic ratio of the ^He nucleus. 

One can write down exchange Hamiltonian for the first and second layers of planar solid ^He in the 
following way: 

pairs triangles rectangles 

where sum in first term is going over all pairs of particles, in second term over all triangles and in third 
term over all rectangles consisting of two triangles. 

The expression of a pair transposition operator Pij has been given by Dirac, 

Pij = l{l + aiaj), (3.4) 

where (jj is the Pauli matrix, acting on the spin at the position number i. The operator P^^ in general 
works in entirely different way, but in case of n = 2 the pair transplonation operators are equal {Pj^^ = Pij)^ 
that we can't write in case of n = 3. 
For n = 3 we have 

-Piife = ^(1 +a-i<7j)(l +CriC7A;), (3.5) 

and 

Using the identity 

((TiC7j)((7iO-fc) = {ajUk) + cri[aj x ak], (3.7) 
we can write the former expression as 

Pijk = ;^(1 + o-iO-j + UjGk + CTkai + ai[aj x a^]) (3.8) 

and 

Pijk = ^(1 + <^i<^j + '^j<^k + CTkCTi + <^iWj X 0-fe]) (3.9) 

hence 

Pijk + P^l = 2 (1 + '^i^j + '^j^k + o-feO-j)- (3.10) 
The four-spin permutation operators can be written as: 

Pijkl = Pijk ■ Pill (3-11) 



Pijkl + (Pijkl) ^ ~ 4 ( + X^"" ' ^f) + Gijkl j , 



where the sum is taken over six distinct pairs (/Ui/) among the four particles (ijkl), and 
So, in terms of Pauh matrices, we have: 



We 



J2 



J3 



{id) 



For the first and second sohd layers, we consider a recursive lattice, instead of the periodic triangular 



(3.12) 



(3.13) 



(3.14) 



one. This lattice is given in Fig. 3.1 




Figure 3.1: The recursive Bethe-type lattice of 4-polygons with additional inner bond. Sq"^ are the spin variables of 
0-th shell, 5*1 of the first shell. 



We can attach to each site of the central plaquette a new one. Carrying out this procedure successively 
for each new shell, we can obtain a recursive lattice which actually is Bethe-type lattice of square plaquettes 
with additional inner links . It is evident that, for each plaquette, the coordination numbers of its sites are 
6 or 5. Introducing the following parameters, 

ai = (3 



"3 = - 

h = (3 



4 ' 



Ja 
4 


4) 


(3.15) 


J4 
4 




(3.16) 






(3.17) 






(3.18) 



the Hamiltonian can be rewritten as: 

-pn = J2[al{S^■Sj + Sj■Sk + Sk■Sl + Sl■Si) + a2{Si■Sk) (3.19) 



+ asiSj ■ Si + {Si ■ Sj) {Sk ■ Si) + {Si ■ 5,) {S, ■ S^) 

- {Si ■ Sk) {Si ■ S,)) 
+ + S| + S| + 5f)}. 

If we use the multisite interaction Ising model, Eq. ( |3.1S| ) takes the foUowing form: 

- PH = {siSj + SjSk + SkSi + siSi) + Q2SjSfc + as {sjSi + SiSjSkSi) 



+h{si + Sj + Sk + si)}, (3.20) 



where s,- takes values ±1 



The partition function takes the form |8C]: 



X9N (4'0 9N (4'^) (4'^) 9N (4'^) , (3.21) 



where Sq"^ are spins of central plaquette, ^Sq°^^ denotes contribution of branch at a-th site of central 
plaquette and N is the number of generations {N — > cxd case corresponds to the thermodynamic limit and 
neglecting the surface effects). For both values (±1) of Sq^\ one can easily calculate: 

gN{+) = a%^Sg%_^{+) + 2bc~^dgl_^{+)gN-i{-) (3.22) 
+ b-^dg'l,_^{+)gN-i{-) + a-%c^d-^gN-i{+)gl_i{-) 
+ 2h-^d-^gN-i{+)gl^i{-) + b-^d'^9N-i{-), 

gN{-) = b-^d^g%_,{+) + 2b-'dgl,_^{+)gN-i{-) (3.23) 
+ a-%c^dg%_^{+)gN-i{-) + b'^d~^gN-i{+)glf^i{-) 
+ 2bc-^d-^gN-i{+)gl^A-) + a%c''d-^g%_^{-), 

where the following notations has been introduced: 

a = expai, 6 = expa2, c = expos, (i = exp/i. (3.24) 

Using Eq. ( 3.22| ) and Eq. ( |3.23| ), one can obtain the recursion relation for variable xn 



giv(+) . 
5iv{-) ■ 



Here 



r( ^ ^ + {2B + l)f,^X^ + (C + 2)f,X + 1 

X^ = /(x^.l), /(^)= ;.3^3 + (C + 2)^2,2 + (25 + l);.x + ^- ^3-2^) 



A = exp(/3(4J3 -2J4 -3J2)) (3.26) 

5 = exp(/3(2J3 - J2)) (3.27) 

C = exp(/5J2) (3.28) 

// = exp(2/i) = exp{p-fhH). (3.29) 
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Figure 3.2: Magnetization plateau for temperature T ~ Q.OAmK (a); bifurcation points and period doubling for 
temperature T — Q.OQmK (b). 



3.2 Magnetic Properties of the Antiferromagnetic Model on the Square 
Lattice 



The recursion relation ( 3.25 ) plays the central role in our further investigations. Because it provides 
all the thermodynamic properties of the system. In particular, analogously, recursion relation for the 
magnetization per site, 

m = ^ f,n ' (3-30) 

can be derived. Let us describe the magnetic behavior of our model in a strong magnetic field. For the 2D 
^He films recent experimental measurements predict the following relations between the exchange energies 
Jn on the regular triangular lattice: 

J3> J2> Ja- (3.31) 

The main features of the resulting magnetic behavior of the system under consideration are caused by the 
interplay between ferromagnetic (J3) and antiferromagnetic (J2 and J4) interactions. Also we use the well- 
known relations between exchange parameters {J2, Js, Ja), which have been estimated from susceptibility 
and specific-heat data in the low density region [|7^ : 

J = J2 - 2J3 -3mK, K = l.STimK. (3.32) 

We have taken the following exchange parameters: J2 = l.lftmK, J3 = 2.35mK and J4 = 1.5mK. At 
T = O.OAmk, we get the magnetization plateau [see Fig. |3]^(a)]. The plateaus at m = and m/irisat = 1/2, 
bifurcation points and period doubling take place for the values J2 = 2mK, J3 = 0.5mK and J4 = 0.5mK 
at T = OMmk [see Fig. 0(b)]. 



3.3 Recursive Approximation to Kagome Lattice 

Since the density ratio of the third layer of "^He is less then the first and second ones, it can use a Kagome 
lattice [p!23|, p.30|]. For using the dynamic system approach it is necessary to approximate the Kagome lattice 



by recursive one. We use the Husimi lattice with two triangles coming out from one site (see Fig. 3. 3). 
In this case 

-He. = J2Y1 (^2 + Yl (^3 + ^3"') (3-33) 

pairs triangles 




and 

nz = -Y^ ^fiH ■ cTi (3.34) 

i 

where 7 is the gyromagnetic ratio of the '^He nucleus. The two-, and three-spin exchanges are given by 
Eq.(3.4) and Eq.( |3.10| ). We have done further approximation passing to classical Ising model as in the ^He 



solid case [B0|. 

The partition function for recursive lattice is written in the form: 

Z = ^ exp(/3/icro) • gl{(To). (3.35) 

Here at takes the values ±1, ctq is the central spin and gn{cro) denotes the contribution of each branch of 
the partition function. 

Introducing, as in previous section, gni+), 5n(— ) and Xn = 9n(+)/5n(— ) one can receive the exact one 
dimensional rational mapping for partition function 

Xn = f[Xn-l), f{x,H,z)= o o i o ' (^-^^^ 

where z = e'^^-^, // = e^t^^ and J = ^-^^2^. 

The magnetization m Eq.(|O0D function of x, temperature and external magnetic field can be 



written as: 

m = . 3.37 

As an example, we take the temperature T = ImK and obtain the figures of m (magnetization) ver- 
sus /i (external magnetic field) (see Fig.( |3.4D ). We take J > since the antiferromagnetic pair exchange 
interaction may be larger than the ferromagnetic three-spin interaction on a Kagome-type lattice. 

At J=0.5mK and 4mK, we have a usual behavior for m (see Fig.( p.4| a) and Fig.( |3.4| b)). The magnetization 
plateau (m=l/3) takes place at J=8mK, Fig.( |3.4| c). At higher values of J, J=18, bifurcation points and 
one-period doubling point occur (see Fig.(^d)) 
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Figure 3.4: Magnetization processes for the values of exchange constants J at temperature T = ImK: (a) J=0.5inK, 
(b) J=4mK, (c)J=8mK and (d) J=18mK 



3.4 Hexagonal Recursive Lattice as an Approximation of the Triangular 
One 

One can write down the exchange Hamiltonian for planar solid '^He in the following way: 

He. = J2Y. + P2^) -h + ^3"') + + P^^) 



pairs triangles rectangles 

J5 {p^ + p^^) + h Y (Pe+p^')^ 



(3.38) 



pentagons hexagons 

where the sum in the first term is going over all pairs of particles, in the second term over all triangles and 



so on (see figure 3.5). 



Using the same technique as in Sec.^j we can derive expressions for fith and sixth exchange ineractions 

P.,fc;m + ^^7fc\„ = ^(l + EK^-)+ Y G^^-V) (3.39) 

t^<'^ fi<u<X<p 
Pijklmn + Pijklmn ~ IQ^^ ^ ^('^/i'^!^) + Yl ^fJ-vXp + Sijl 



'ijklmn ) 



(3.40) 



p.<u<X<p 



where Sijkimn is 

Sijklmn 



= [{criaj){akCri){(7mO-n) + ('7jCrfc)(cr/(7m)((T„crj)] 

- [(0"i(Tj)((Tfc0-m)((T;(T„) + (cTjCTfc ) ((T/(7„) (fTmlTi ) 

+ io-kCri){amCri){crnaj) + (cri(Jm)((T„crj)((Ti(Tfc) + (crmCr„) 

+ (o"i(Tfc)((TjCrz) + (cr„crj)((Tj((T„crfc)cTm)] + [(o'j(Ti)(crjCr„)(crfcCrm) 

+ ('7j(Tm)(<7fcfTj)((T/(T„) + ((Tfc(T„) ((T/CTj ) ((Tm(Tj )] + [(cTjCrj ) ((Tfc(J„) ((T/CTm) 

+ {o'jCrk){criai){amO-n) + (<7fc(T;)((TmCrj)(cr„(Tj)] 

- [{(^i(^l){(^j(^m){(^k(^n)] 



(3.41) 




Figure 3.5: The hexagon recursive lattice. 



Now let us consider the recursive lattices which are connected through the sites for getting one— dimensional 
rational mapping which is close to triangular lattice. The central hexagon is approached to the triangular 
lattice by connection of center to vertices of hexagon, so we get six triangles in the central hexagon. The 
recursive lattice is constructed by adding next shells (branches) of hexagons on vertices of previous one (see 
figure |3.5D . Next shells are constructed in the same way. For hexagonal recursive lattice we obtain the exact 
recursive relations for partition function. It's important to mention that this recursive relation is dynamic. 

By now we considered Heisenberg model for description of our recursive lattice. We can assume that 
in the strong external magnetic field all atom' spins will be directed mostly by certain z axis. So it 
will be possible to calculate the eigen-values of spin— operators. One may consider the variables Si as 
classical vectors as well as Pauli matrices. If we use the multisite interaction Ising model, we can take 
+1 and -1 instead of cjj. In the rest of the paper all Constance will be taken in temperature's scaling 
(J2 =J2/r, Js =J3/T, J4 =J4/r, Js =J5/T, =Je/T, h =H/T,). 

When the lattice cut apart from the central hexagon, it separates into six identical branches. So we can 
first realize a summation over all spin configurations on each branch, getting the same result every time, and 
then the sum over spins of the central hexagon. Here follows, that we can represent the partition function 
as 

p~h{so+Si+S2+S3+S4+S5+Se)+Hex 



X9n{si)gn{s2)9nis3)gn{s4)9n{s5)gn{s6) 



(3.42) 
(3.43) 



from which the expression takes the form 



where 



Z = e- 
+ 6e 



+ 



+ 6e-3'^-«-^2+6-^3-4J4+2J-5^(+)4^^_)2 ^ 2e-'^-3-^^5(+)35(-)3 

+ 6e-'*-^-^2+^-^3-2J4^(+)3^^_)3 ^ 6e'^-^-^2+4J3-2J4^(+)3^(_)3 
+ 9e'^-4^2^(+)2^^_)4 ^ 6e'^-6-^2+2J3^(+)2^^_)4 ^ 3g3/.-6J2+4J3^(+)2^^_)4 
+ 6e3'*-6-^2+^-^3-2J4^(+)2^^_)4 ^ gg3^-8J2+6J3-4J4+2J5^(+)2^^_)4 
+ 6e=^'*-5'^2^(+)3(-)5 + 6e^''-''-^2+8J3-6J4+4J5^(^^)^(_^5 ^ g5ft-6J2-2J6^(_)6 
7/1-12 J2+I2J3-I2J4+I2J6-2J6 , 



y ^ g-fe(si+Sl+S2+S3+S4+S6+S6)+-ffex 

/ (2)\ / (2)n / (2)n / (2)x , (2)x 

XC^n-l(Sl ')5n-l(S2 )5n-l(S3 )9n-l{S4 )9n-l[Si, ')■ 



We can obtain the recurrent relation for the variable Xr 



(3.44) 



(3.45) 



fix) 



g-4/i+7J2 _|_ g-6/i+3J2+8J3-6J"4+4J5 _|_ (^3g-2/i+8J2 



^g-4W2+4J3 ^ 2e 
+ (e9-^2 _^ g-2/i+9J2 _^ gg7J2+2J3 _^ gg-2/i+7J2+2J3 
_|_3g5,/2+4J3-2J4 _|_ 3g-2?i+5./2+4J3-2J4-)j^2 

+ (^Qe^-h + 4g6J2+2J3 _^ 2e2^+'^"'2+4J3 

_|_4g2/i+6J2+4J3-2J4 _(_ 4g2/i+4J2+6J3-4J4+2Js^ v3 
+ (5e^'*+^'^2 _|_ 5g4/i+3J2+8J3-6J4+4JB j^4 



-4/1+6J2+4J3-2J4 



26-2^1+6^2+273 
2g-4/i+4 J2+6 J3-4J4+2 Jb ) 



')X3 



(3.46) 



_|_ ^g4/i+6J2-2J6 _|_ g6/i+12J3-12J4+12J5-2J6^j^5 



^g4/i+7J2 _|. g6/i+3J2+8J3-6J4+4J5 



)X^ 



+ (3e^'^"^^"^2 + 2e^'*+^"^2+2^3 _|_ g4fe+6J2+4J3 

_|_2g4/i+6J2+4J3-2J4 _|_ 2g4/i+4J2 +673-474+276-1^^4 
^3g5J2+4J3-2J4 _|_ 3g2fe+5J2+4J3-2J4-)j^3 

+ (ee^'^^ _^ 4g6J2+2J3 ^ 2e-2'^+^-^2+4J3 _^ 4g-2/l+6J2+4J3-2J4 

_^4g-2fe+4J2+6J3-4J4+2J5-|j^2 _^ |^gg-2/i+7J2 _^ 5g-4;i+3J2+8J3-6J"4+4J6-)j^ 
+ g-4/i+6J2-2J6 ^ g-6/i+12J3-12J4+12J5-2J6 



This recurrent relation plays a crucial role in our further investigations, because it's getting possible 
obtain all thermodynamic quantities of the considered system. 



3.5 Magnetic Properties of the Antiferromagnetic Model on the Hexag- 
onal Recursive Lattice 

Using the technique of dynamic system theory, apphed to one— dimensional rational mapping, we can 
get magnetization relations per site. As mentioned above for an homogeneous lattice the magnetization 
function is given by the formula 



E 



m = 



is) 



-(3H 



(3.47) 



E(.) e- 

where Si is the spin on an arbitrary site. In our case, the lattice is homogeneous, so on central hexagon we 
divided the tree of hexagons on six branches. 

We could obtain the magnetization value of sublattice composed by hexagon central vertices. 

1 



z E 



and for magnetization sublattice composed by corner verticies. 
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- E 
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So from this expressions finally we could get 
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-5fe-6J2-2J6 



- 7/1-12 J2+I2J3-I2J4+I2J5-2J6 



(e ""^ — e 

gg-3/i-5J2 _ gg-5/i-9J2+8J3-6J4+4J5)j^ 

gg-fe-4J2 _^ gg-/l-6J2+2J3 _ 3g-3/l-6J2+4J3 

_gg-3/j-6J2+4J3-2J4 _ gg-3/i-8J2+6J3-4J4+2J5)jf2 
2gft.-3J2 _ '^^-h-ZJi _|_ Y2^^-^Ji+2Jz _ -[^2g-''--5J2+2J3 

_^gg/l-7J2+4J3-2J4 _ gg-/l-7J2+4J3-2J4)J^3 
gg/i-4J2 _^ gg?i-6J2+2J3 _^ 3g3/i-6J2+4J3 

_^gg3/i-6J2+4J3-2J4 _|_ gg3/i-8J2+6J3-4J4+2J6)j^4 
gg3/i-5J2 _|_ gg5?i-9J2+8J3-6J4+4J6)j^5 

-7/1-12 J2+I2J3-I2J4+I2JB-2J6) 



(3.50) 



g5?i-6J2-2J6 _|_ g7/i-12J2+12J3-12J4+12JB-2J6^ 



-5/1-6 J2 -2 Je 



+ e 



gg-3?i-5J2 _^ gg-5/i-9J2+8J3-6J4+4JB)j^ 



9e 

+6e- 



-/1-4J2 



+ 6e" 



-/1-6J2+2J3 



+ 3e 



-3?i-6J2+4J3 



+ 6e 



-3h-6J2+4J3-2J4 



-3/1-8J2+6J3-4J4+2J5- 



2g-ft.-3J2 _j_ 2g't-3J2 _j_ -|^2e^''^5'^2+2J3 _|_ -[^2g'*~5J2+2J3 



+6e- 



-fc-7J2+4J3-2J4 



+ 6e 



h-7J2+AJ:i-2J4^ 



gg/i-4J2 _^ gg/i~6J2+2J3 _^ 3g3fe-6J2+4J3 
_|_gg3h-6J2+4J3-2J4 _|_ gg3/i-8J2+6J3-4J4+2J5)j^4 
gg3/i-5J2 _|_ gg5/i-9J2+8J3-6J4+4J6)j^5 
g5h-6J2-2J6 _|_ g7/i-12J2+12J3-12J4+12JB-2J6)jf6j 



and 

j^^ ^ |-_g-5/i-6J2-2J6 _ g-7/i-12J2+12J3-12J4+12J5-2J6 ^3 ^^^^ 

- (4e~^''~^'^2 _^ 4g-5/i-9J2+8J3-6J4+4J5>|Jj^ 

- (Se'^"'^-^^ + 2e~''~^'^2+2J3 _^ g-3/i-6J2+4J3 

_|_2g-3/i-6J2+4J3-2J4 _|_ 2g-3/i-8J2+6J3-4J4+2J5-jj|^2 

+ (3e''"*^'^2 _^ 2e^"^'^2+2J3 _^ g3?t-6J2+4J3 

_|_2g3/i-6J2+4J3-2J4 _|_ 2g3/i-8J2+6J3-4J4+2J5-jj^4 
+ (-4g3/i-5J2 _^ 4g5/i-9J2+8J3-6J4+4J5~|j|^5 
_^ |^g5/i-6J2-2J6 _^ g7/i-12J2+12J3-12J4+12J5-2J6^j^6j 

|-^g-5h-6J2-2J6 _|_ g-7/i-12J2+12J3-12J4+12J5-2J6-j 

_|_|'gg-3/i-5J2 _|_ gg-5h-9J2+8J3-6J4+4J5-jj^ 

_|_ (9g-'^-4-'2 _|_ gg-/l-6J2+2J3 _(_ 3g-3/l-6J2+4J3 

_|_gg-3/i-6J2+4J3-2J4 _|_ gg-3/i-8J2+6J3-4J4+2J5-jj|^2 

_|_gg-'l-7J2+4J3-2J4 _|_ gg/l-7J2+4J3-2J4-)^3 

+ (9e^~^'^2 + 6e^~^'^2+^'^^ + 3e3''-6-^2+4J3 

_^gg3h-6J2+4J3-2J4 _^ gg3?t-8J2+6J3-4J4+2J5^jl^4 
_|_ |-gg3/i-5J2 _|_ gg5/i-9J2+8J3-6J4+4J5-jj|^5 
_|_ |^g5/i-6J2-2J6 _|_ g7/i-12J2+12J3-12J4+12J5-2J6^j^6j 

Having these dynamical expressions one can draw the plots of magnetization vs. external magnetic field 
for sublattices in various temperatures. To do that one should fix the value of the dimensionless magnetic 
field h (for a given temperatures and exchange parameters) and implement the simple iteration from the 
recursion relation for /(x), beginning with some initial xq. For achieving to the thermodynamical limit we 
have to apply infinite amount of iterations (n — > oo). 

From experimental measurements and theoretical calculations we could conclude that the relations be- 



tween the Jn exchange energies on the regular triangular lattice are||68||[70|[69|: 

J3 > J2 > J4 > Je > J5 (3.52) 

It's important to mention that the values of pure J„ are not observable in the experiment measurements, 
because each n-spin exchanges makes also a contribution to a few (n — l)-spin exchanges, and thus there 
are some effective exchanges parameters, which are certain contributions of J„ (such as J = J2 — 2J3 and 
K = — 2J5) and can be directly obtained from experiments [|^, 70, |128| . 

The exchange parameters depend on the particle density and on the type of the lattice, particularly 
on its coordination number and dimensionality. For instance, it was found that for 2D triangular lattice 
at high densities the exchange J is dominant, mainly due to the three-spin exchange, but the ratio \K/J\ 
increases rapidly with the lowering of the particle density, and below other exchanges become important 
also. The magnetic properties of the system (ferromagnetic or antiferromagnetic) depend on which exchange 
interactions are dominant at the present value of the particle density. 

From the arguments stated above one can conclude that there is a large freedom in the choice of concrete 
values of the exchange parameters J2, J3, J4, J5 and Jg. Moreover, it is quite difficult to identify our model, 
having so many assumptions, with some concrete value of particle density. 
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Figure 3.6: J3 = 2.5, J5 = 0.5, J2 — Ja = Jq = Q pure ferromagnetic case, a) is the magnetization for the central 
vertex lattice with temperature T = 0.1, b) is the magnetization for the corner vertex lattice with temperature T — 0.1, 
c) is the average magnetization for the lattice with temperature T = 0.1 



The simplest case of pure ferromagnetic behavior ( J3 = 2.5, J5 = 0.5, J2 = J4 = Je = 0) is presented on 
figure Figure 3^ a represents the magnetization for the central vertex of hexagon recursive lattice with 
temperature T =0.1, figure 3.6b the magnetization for the corner vertex lattice with temperature T = 0.1 



and figure |3.6| c the average magnetization for the lattice with temperature T = 0.1. At relatively high 
temperatures the magnetization curve has a smooth monotone form of Langevin type with a rather large 
value of the saturation field. Decreasing the temperature, the curve becomes more steep and the value of the 
saturation field decreases. Further decreasing the temperature leads to the magnetization diagram. We can 
see smoothly increasing curves, which are expected and corresponding to only ferromagnetic interactions. 

The physical case (J = 3mK) with ferromagnetic and antiferromagnetic interactions ( J2 = 2, J3 = 
2.5, J4 = 1.8, J5 = 0.5, Je = 1) for enough high temperature T = 10 is represented on figure |3.7| . On plots 
we can see the paramagnetic limit as expected in high temperatures. 

There is represented the depends between external magnetic field and magnetization on figure (a) 




Figure 3.7: J2 = 2, J3 = 2.5, J4 = 1.8, J5 = 0.5, Je ~ 1 paramagnetic limit at enough high temperature, a) is the 
magnetization for the central vertex lattice with temperature T = 10, b) is the magnetization for the corner vertex 
lattice with temperature T ~ 10, c) is the average magnetization for the lattice with temperature T = 10. 



is the magnetization for the central vertex lattice with temperature T = 0.1, (b) is the magnetization for 
the corner vertex lattice with temperature T = 0.1, (c) is the average magnetization for the lattice with 
temperature T = 0.1 for physical case with J = 3mK (J2 = 2, J3 = 2.5, J4 = 1.8, J5 = 0.5, Jg = 1). The 
hatch line represents the magnetization value with m/msat = 2/3. On figure we can see bifurcation 
points in high magnetic fields. This fact is making clear the advantage of represented model in high external 
magnetic field. There is a one period doubling behavior between bifurcation points, so that the contribution 
of antiferromagnetic interaction is visible. 

The magnetization plateau on m/ ragat = 2/3 on figure |3.9| c (a) is the magnetization for the central vertex 
lattice with temperature T = 0.1, (b) is the magnetization for the corner vertex lattice with temperature 
T = 0.1, (c) is the average magnetization for the lattice with temperature T = 0.1) with (J2 = 2, J3 = 
2.05, J4 = 0, J5 = 0, Jg = 0). As in previous example we can see bifurcation points in high magnetic fields 
and a one period doubling behavior between bifurcation points. 



From figures |3.9| a, |3.9| b and |3.9| c we can conclude that in represented recursive lattice it's possible to 
separate two possible kinds of lattice's structure with four sublattices (shells of the considered hexagon tree) 
each one represented on figure gl0| and |1| which are changing with each other the directions of their's 
spins positions during the process of interaction. That phenomena is called modulated phases process. 
The Roman numbering is corresponding to each sublattice and the bottom side of each figure represents 
one period of structure. We can conclude from doubling of plots on figures ^ a and 3^b the Heisenberg 
interaction between spins of hexagonal recursive lattice are going on opposite phases depending on the fact 
of dominating either ferromagnetic or antiferromagnetic essence of interaction. 

Supporting on results reported above we are able to conclude that the whole lattice maybe constructed in 



two ways of recursively repeating shells represented on bottom side of figure |3.10| and figure |3.11| . The possible 
not monosemantic behavior of lattice structure in the thermodynamic limit may come from difference of 
initial value of ratio xq = g{+)/ g{—)- 

In the following plots (J2 = 2, J3 = 2.3, J4 = 1, J5 = 0.1, Jg = 0.5) we can see bifurcation points and one 
period doubling phenomena between bifurcation points represented on figure |3l^ (a) is the magnetization 
for the central vertex lattice with temperature T = 0.8, b) is the magnetization for the corner vertex lattice 
with temperature T = 0.8, c) is the average magnetization for the lattice with temperature T = 0.8). 

Thus we've got magnetization plateaus, bifurcation points, doubling behavior, modulated phases for 
multiple exchange model where is taken into account both ferromagnetic and antiferromagnetic interactions. 
When we are changing interaction parameters or as the same the density in the physical limits got from 
experiment [70|] (close to J = 2>mK) we see magnetization plateaus, bifurcation points, doubling behavior, 
modulated phases at low temperatures, (fig. 3^,3.£. 3!T^ ). 




Figure 3.8: J2 = 2, J3 = 2.5, J4 — 1.8, J5 = 0.5, Jg = 1 physical J = ZmK case with ferromagnetic and 
antiferromagnetic interactions contributions. There are bifurcation points and one period doubling between bifurcation 
points, a) is the magnetization for the central vertex lattice with temperature T = 0.1, b) is the magnetization for 
the corner vertex lattice with temperature T = 0.1, c) is the average magnetization for the lattice with temperature 
T = 0.1. 




Figure 3.9: J2 = 2, J3 = 2.05, J4 = 0, J5 = 0, Jg = case with ferromagnetic and antiferromagnetic interactions 
contributions. Tlicre are bifurcation points, a one period doubling between bifurcation points and magnetization 
plateau at m/msat = 2/3. a) is the magnetization for the central vertex lattice with temperature T = 0.1, b) is the 
magnetization for the corner vertex lattice with temperature T = 0.1, c) is the average magnetization for the lattice 
with temperature T = 0.1. 
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Figure 3.10: First of the two possible kinds of lattice's structure with four sublattices (shells of the considered hexagon 
tree). 




Figure 3.11: Second of the two possible kinds of lattice's structure with four sublattices (shells of the considered 
hexagon tree). 




Figure 3.12: J2 = 2, J3 = 2.3, J4 = 1, J5 = 0.1, Je = 0.5 physical case with ferromagnetic and antiferromagnetic 
interactions contributions. There are bifurcation points and one period doubling between bifurcation points, a) is the 
magnetization for the central vertex lattice with temperature T — 0.8, b) is the magnetization for the corner vertex 
lattice with temperature T — 0.8, c) is the average magnetization for the lattice with temperature T — 0.8. 



Chapter 4 

Face-cubic Model i 



Cubic symmetry plays a significant role in many types of phase transitions and critical phenomena. 
Still in 70-s various models possessing cubic symmetry have been introduced to determine the nature of the 
displacive phase transitions in perovskites |131, [132| . Magnetic phenomena in cubic crystals are also affected 
by the lattice structure. For instance, in the crystalline solids with cubic-symmetric lattices (Fe, Ni, etc.) it 
leads to the modification of the magnetic exchange interaction giving rise to additional contributions to the 
conventional 0(N)-symmetric Heisenberg Hamiltonian. The simplest contribution of the underlying cubic 
symmetry is the single- ion anisotropy of the form X^jS^. During the last three decades different aspects of 
this issue have been the subject of various investigations. 

In the framework of the field-theoretical approach to critical phenomena modelling of effects of cubic 
symmetry is usually performed in terms of continuous-spin Landau-Ginzburg effective Hamiltonian with 
cubic anisotropy, i.e. the c/)'^ - theory with an additional cubic term which breaks explicitly the 0(N) 
invariance to a residual discrete cubic symmetry [133| - | 138 |: 



+ -a^' + -I 



N 



(4.1) 



Q = l 



where (^{x) = (01, (^2, ■ ■ ■ , <Pn) is a continuous N-component local vector order parameter, 0^ = (0^)^ and 
f3 = 1/kT is inverse temperature. The spin orientations, orthogonal to faces of an A^-dimensional hypercube 
will be favorable for v < 0, whereas positive values of v favor the orientations toward the corners. 

This model received numerous important applications, among which are such as the oxygen ordering in 
YBaCusOg+x [134], one of the most studied high-Tc superconductor; the buckling instability of confined col- 



loid crystal layer |135U ; the micellar binary solution of water and amphiphile |136| etc. The two-dimensional 
case has been recently intensively studied within the framework of renormalization group technique in the 



space of fixed dimensionality up to five-loop approximation [137, 138]. The model was found to have four 
fixed points: the Gaussian one, the Ising one with N decoupled components, the 0(N)-symmetric and the 
cubic fixed point. Analogous calculations were made for D = 3, as well revealing some peculiar type of cubic 
fixed point for N > 2 corresponding to specific anisotropic mode of the critical behavior [139|-[142]. The 
surface critical behavior of the model was also studied ]14S]. It has been recently shown that the behavior 
of the spherical many-spin magnetic nanoparticle with surface anisotropy may be modelled by an effective 



single macro-spin with cubic anisotropy terms in the effective energy ]144]. 

Another class of lattice cubic-symmetric spin models was introduced by Kim, Levy and Uffer in Ref. 



]145] to explain the tricritical-like behavior of cubic rare-earth compounds, particularly holmium anti- 
^The results considered in this chapter are published in Ref. [Bl[. 



monide, HoSb |145]-||T5^]. Projecting the pair exchange interaction onto the sixfold degenerated ground-state 
manifold of Ho^^ ion they arrived at the following effective Hamiltonian: 



I3n 



(4.2) 



where the classical spin variables Sj are the unit vectors restricted to have orientations, orthogonal to the 
faces of the cube and the sum is going over all the pairs of nearest-neighbor sites. Generalization for Q 
component spin is obvious: each spin can assume 2Q orientations: 



S, G {(±1,0,...,0),(0,±1,...,0)...(0,...,0,±1)}. 



(4.3) 



This model is known as the Face-cubic model. It is connected to the continuous cubic model of Eq. (5T) via 
the limit of strong anisotropy {\v\ ^ |5|). In a similar way one can also consider quadrupolar pair interaction 
terms and obtain the following Hamiltonian provided all interactions are homogenous [Il4'i 



141 : 



/3Wfc = j5^(SiS,) + ir^(S,S,)2. 



(4.4) 



It is easy to see that Hamiltonian ( |4.4| ) can be represented in terms of two sets of discrete variables: a 
Potts-like determining which component of Sj is non-zero and an Ising-like, corresponding to the sign of 
the component. Indeed, 



where (Tj = ±1 and ai = 1,2, ...Q. Therefore we have 



{hi) {hi) 



(4.5) 



(4.6) 



Formally we can enlarge the Hamiltonian ( [4.4| ) up to the interaction terms of higher power of the (SjSj): 



{i,i) ^=0 



(4.7) 



For arbitrary finite L this expression leads to the same Hamiltonian that of Eq. (4.6) with J 
and K = X]fc=o'^2fc- At = the Hamiltonian ( ^^ ) reduced to the 2Q - state Potts model; J = 
corresponds to two decoupled Q - state Potts model, and at Q = 2 one obtains the Ashkin-Teller model. 
Variational renormalization-group study of the pure and diluted Q - component face-cubic model in two 
dimensions has revealed existence of four competing possible types of critical behavior corresponding to Q 
- state Potts model, 2Q - state Potts model, Ising model and special "cubic" fixed point |14£, |15C[| . It was 
also found that at Q < Qc = 2 the transitions are continuous and critical behavior of the discrete face-cubic 
model belongs to the 0(N)-model universality class. For Q = 2 the Ashkin-Teller-like behavior occurs and 
for Q > 2 transitions are of first order. There are also early results obtained within the mean-field (MF) 



theory |145| ] and using the Bethte-Peirels (BP) approximations and high-temperature series |146| for the 
case when only dipolar pair interactions are included {K = 0). In the BP approximation the critical value 
of spin component Qc, above which transitions are of the first order was found to be given by 



Qc 



2 

1 + 3, 



1 + 5 

q 



1/2 



1 



(4.8) 



Figure 4.1: The Cayley tree with coordination number q = 3 and 3 shehs. 



where q is the coordination number of the lattice. The Umit q ^ oo corresponds to the MF solution and 



gives Qc = 3. In Ref. |14C] the high-temperature series for Q ~ component face-cubic model on three 
dimensional fee lattice were constructed up to 5-th order, from which authors obtained Qc = 2.35 ± 0.2, 
whereas Eq. (|4.8|) gives us Qc — 2.8 Moreover, although the MF solution of cubic model predicted the 



tricritical like behavior |145|] , the BP approximation showed that it is not the case |146|. Only inclusion of 



single-ion-anisotropy terms, quadrupolar pair interactions and crystal fields may drive the system tricritical 



4.1 Recursive Methods, Cayley Tree and Bethe Lattice 

Among the vast variety of statistical mechanics lattice models with strong local interactions only very 
limited amount allows exact solutions. These exact solutions are known only for low-dimensional systems, 
more precisely for d = 1 and d = 2; most of exact solutions for two-dimensional systems are known only 



in the absence of the external field and/or for the special choice of model parameters [92|. Well known 
conventional approximate methods like MF theory and BP approximation in general can provide only more 
or less qualitatively satisfactory picture and in some cases they just fail. That is why the quest for the 
alternative approaches, which can provide more reliajale results for the thermodynamics of lattice models is 
very important. 

The recursive lattices are of twofold interest. On the one hand, the models, defined on a recursive lattices 
can be considered as an independent original area of research in which the powerful methods of dynam- 
ical systems theory and fractal geometry are successfully exploited for determining the thermodynamical 
properties of the statistical mechanics models. On the other hand, recursive lattices provide a specific kind 
of approximate treatment of many-particle systems physics. In the heart of these approximations not the 



simplification of tlie character or /and strength of the interaction between the system elements, but the 
modification of the topology of the underlying lattice lies. This modification in most cases consist in the re- 
placement of the regular periodic "physical" lattice by a recursive one, constructed by the certain algorithm 
and possessing the self-similarity, maintaining all interaction unchanged. For the large class of lattice mod- 
els with commuting variables ("classical models") this approach leads to the exact solution of the statical 
mechanics problem in terms of the theory of dynamical systems or, more precisely, in terms of discrete maps. 
The solution formally is rather similar to that which corresponds to conventional BP approximation but, in 
contrast to the latter, it is the exact solution for the lattice endowed with the recursive structure provided 
the boundary sites are properly taken into consideration. It was argued that in some cases, particularly in 
the models where multi-site interactions are presented 87, ^] the recursive lattice approximation gives 
more reliable results than conventional BP [B6|. 

The simplest example of recursive lattices is a Cayley tree, the self-similar graph, which contains no 
cycles. In order to construct a Cayley tree one should start with a central site O. At the first step this 
site should be connected by links with q others, which constitute the first shell or the first generation of 
the Cayley tree. The second shell of the Cayley tree is constructed by repeating this procedure for all sites 
of the first shell: each of the q sites belonging to the first shell connects to the q — 1 new sites. Thus, the 
second shell contains = <l{(l — 1) sites. Accomplishing this construction for n steps one arrives at the 
recursive connected graph, which contains no cycles and is called the Cayley tree with coordination number 
q and n generations (shells) (Fig. |4.1| ). As the number of sites in A:-th shell is q{q — 1)^^^ the total number 
of sites then 

iVn = 1 + E g (g - 1)^-^ = (4.9) 

Thus, each site of the Cayley tree has coordination number equal to q except the sites on the last shell 
which have only one neighbor. The peculiar property of the Caylet tree is the large amount of boundary 
sites which even in the thermodynamical limit n — > oo comprise a finite fraction of the total number of sites 

q{q-lf (9-2)g(g-l)"-^ q-2 

^™ AT = ^™ 1 5 = 7- 

n^oo JSn n^oo q [q — 1) —2, 9 — 1 

This feature causes anomalous properties of the statistical mechanics systems on the full Cayley tree. For 
instance, in Ising model there is no zero field magnetization, whereas the derivations of free energy with 



respect to the external field exhibit singular behavior |151] - [154|. 

The Bethe lattice is an object intimately linked to the Cayley tree but, in contrast to the latter, it is 
devoid of the complications which arise from the boundary sites. When one deals with the Bethe lattice the 
general structure and topology of the Cayley tree are preserved whereas only contributions from the bulk 
sites lying deep inside the Cayley tree are taken into consideration. In this regard the Bethe lattice is the 



interior of the Cayley tree |9^, 155]. Dealing with the Bethe lattice we suppose the underlying Cayley tree 
to be large enough to achieve thermodynamical limit and consider only the sites situated far away from a 
boundary. One can regard these sites as a uniform lattice with coordination number q. 

An undeniable advantage of the Bethe lattice is the possibility of exact solutions for many types of 
statistical mechanics problems. A large amount of problems of the theory of magnetism, macromolecule 
physics, lattice gauge theory, self-organized criticality, dynamical mean-field theory (DMFT) and general 
questions of statistical mechanics like zeroes of partition functions have been successfully considered on 
Bethe or Bethe-like lattices, revealing many interesting exact results and deep connections between the 
theory of dynamical systems and statistical mechanics [151] - ]164]. 



4.2 Face-cubic Model on Planar Graphs 



It is possible to calculate exactly the partition function of the FCg-model on Cayley tree in case of 
the absence of external field. Let us consider the partition function of the model without external field on 
arbitrary planar graph G 

Zg = Y.Y. n exp{(Ja,c7,+if)5«,,«J. (4.11) 

{a} {a} (i,j)&G 

As in the case of Potts model [|165|| one can represent it in the following way: 



^G=EE n (i+f^<x../a.,a,), (4.12) 

where L^o-io-j = e'^"'-''^^^ — 1. Then summing out over the all ai we obtain that 

^G=Y.Y.Q'''^'''^ n f^-^.' (4-13) 

{a} G'CG {i,j)&G' 

where the second sum is going over all spanning graphs G' of the underlying planar graph G and Rq{G') 
is the 0-th Betti's number of G' which coincides with the number of connected components of G' and the 
product is over all links belonging to G' . Using the identity aiaj = 2(5o-.o- — 1 one can represent Eq. 



(4.13) in terms of double power series associated with the underlying lattice provided latter is a planar 



graph: 



G'CG {a} {i,j}eG' 

u = e^-'-l, 
e'^J - 1 

(4.14) 



where e(G') is the edges number in the graph G' and all cjj assume values or 1. So, for each spanning 
subgraph G' of planar graph G we obtain a partition function of 2-state Potts model: 



G'CG 

= Yl Q'^'^'^'^u^'''^ Yl 2««(^'%^(«"), (4.15) 
G'CG G"CG' 

where the second sum is going over all subgraph of G' Thus, we have succeeded in representing the partition 
function of the FGq - model given by Hamiltonian ( [4.6D in terms of double power series associated with 
the underlying lattice provided the latter is a planar graph. As in case of ordinary spin Potts model one 



can make l/Q-expansions for FCg-model with the aid of Eq. ( 4.15 ). For instance, for the two-dimensional 
square lattice one obtains 



= {2Qf{l + 2N^q + l) + N{2N-l){^fq + lf (4.16) 

+ N^^iil + 1)' + il)') + (CIn - Ni2N - m^fil + If 
+ Ar(2Ar_4)^((^ + l)5 + (|)4 + (^)5) 



Now we can use the well known relations between the topological invariants of graphs |166| ] to rewrite 
Eq. (4.15) in another way suitable for our further purposes. Applying the Euler theorem to planar graph 
G we get 



Ro{G)-Ri{G) = n-e{G), (4.17) 

where the first Betti's number Ri{G) coincides with the number of independent cycles of G and n is the 
number of sites. Therefore 

ZG = {2Qr Q^^(^')(^)^(^') 2«^(«")(|)-(G"). (4.18) 

G'CG ^ G"CG' 



Partition function for FCg-model in the form of Eq. ( 4.18 ) can easily be calculated in case of the 



so-called forests, the planar graphs, which contain no cycles. Connected part of forest is called tree. Let us 



suppose that we have a tree T with n sites and L„ edges. Then Eq.( |4.1q ) takes the form 

2t = {2Qr E E (Ir'''"' 

G'CT ^ G"CG' 

Ln k 

= (2Q)"EGt(^)'E^i(i)^ (4-19) 

A:=0 ^ 1=0 

where = are the binomial coefficients and the sum is going over all subgraphs with given volume 

k. With the aid of Newton's binom we obtain 

If T is the Cayley tree then Ln = n — 1 and we will have 

Zcayley = 2Qiuiv + 2) + 2Q)''-\ (4.21) 

Thus, the free energy per site for FCg-model on Cayley tree is 

/ = -ksT lim ^Cayley ^ + 2) + 2Q) 

= - A:bT log 2 - /cbT log (e^ cosh J + Q - I) . (4.22) 
It is noteworthy that this result is in full agreement with that obtained by Aharony for one-dimensional 



model by transfer-matrix technique [148|. Indeed, a one-dimensional chain can be regarded as a "tree" in the 
sense mentioned above. Thus, the free energy of the FCq model on a Cayley tree in a thermodynamic limit 
is continuous for all T. Formally it coincides with that for one-dimensional systems where the continuity is 
really the case. But, as is known, the BP approximations for the lattice model with an arbitrary number of 
component solves exactly the problem on a Cayley tree [|159[|. Therefore the model on a Cayley tree must 



possess a finite critical temperature as predicted by BP approximation [146|. The origin of this peculiar 



properties of the Cayley tree was revealed by Eggarter |151] on the example of Ising model. He argued that 
the equivalence of all sites in the thermodynamic limit, which is one of the key points of BP approximation, 
breaks down for the sites of the Cayley tree which are situated close to the surface. As appears from Eq. 



( CT) these sites comprise rather large fraction of the sites of the Cayley tree and, thus, they determine the 



behavior of all thermodynamic quantities to a great extent. 

However, even though no phase transitions in the thermodynamic limit occur at the value of critical 
temperature Tc predicted by the BP method the value of the order parameter (magnetization etc.) in the 



interior of the lattice, i. e. for the region far from the surface, which is generally called "Bethe lattice", 
will undergo jumps from zero to its BP value when the temperature is below Tc- In order to illustrate it we 
will use the technique of dynamical system theory, which becomes a powerful tool for investigating various 
physical problems on recursive lattices. 

4.3 Recursive Method for Face Cubic model on Bethe lattice. 

Many statistical systems defined on the recursive lattices are famous for the possibility of exact solution 
in terms of dynamical system theory. In the heart of these exact solutions the self-similarity of recursive 
lattices lies. For instance, if we cut apart the Cayley tree at the central site it will give q identical branches 
each of which is the same Cayley tree with the number of generation decreased by 1. Using this fact one can 
establish the connection between the partition function of a model defined on the Cayley tree containing 
n shells with the partition function of the same model defined on the Cayley tree containing n-1 shells. 
Therefore, the thermodynamical problem is reformulated in terms of discrete maps, given by recurrent 
relations. Let us introduce a symmetry breaking field into Hamiltonian ( [4.61 ). The field can be considered 
as a uniform magnetic field pointing along the first coordinate axis H = 0, ...0). The corresponding 
interaction Hamiltonian is of the conventional Zeeman type: 

Wz = -H^Si, (4.23) 

i 

which leads to the following form of the FCg-model Hamiltonian in the field: 

- (3rLFC = (^i(^jSa,,a, +KY^ 6a,,a, + X] (4-24) 

where h = f3H. Hereafter we pass from the Cayley tree to the Bethe lattice having in mind one of the ways 
mentioned in Section II. Thus, we should no longer care about boundary sites and boundary conditions. 
According to that one can represent the partition function of the FCg-model on the Bethe lattice in the 
following form: 

Z= e'^''"'"'"M5n(^o,ao)]^ (4.25) 

(o"0.ao) 

where gq and ao are the variables of the spin in the central site and gnic^o, cto) refers to a partition function 
of individual branch: 

5„(fTo,ao) = Y Y exp(^(Jcro(Ti + K)5ao,ai + ^(JfTiO-j + K)6a,,aj + • 

Each branch, in its turn, can be cut at the site which was previously connected to the central one. 
This will give us q-1 identical branches being the Bethe lattices with n-1 generations. Thus, the connection 
between g-n and gn~i is 

fl'n(o-o,ao)= Yj exp ((JcJoo-i i^) 5 QQ ,Q;i 

)[gn-i{ai,ai)f-\ (4.26) 

Here we obtain the system of 2Q recursion relations but, in fact, only three of them are independent as all 
gn{cF, a) for Q 7^ 1 are identical. So, from (|4.26| ) we have 

5„(+, 1) = e-^+^+''[g„_i(+, 1)]^-^ + e-^+^-"[5„-i(-, l)]"^'^ + 2(Q - l)bn-i(±, (4.27) 
5„(-,l)=e"^+^'+'^[5„-i(+,l)ri + e-^+^-'^[5„-i(-,l)ri + 2(Q-l)b„-i(±,*)r\ (4.28) 
ff„(±, *) = e^bn-i(+, 1)]"-' + e-^b„-i(-, 1)]^-^ + (e^+^ + e-^+^ + 2{Q - 2))[5„-i(±, (4.29) 



where gn{^, *) stand for any (2Q — 2) partition functions corresponding to the individual branch with Si 
whose direction is non cohinear with the first coordinate axes. Introducing the variables 

<?n( + ,l) 



5n(±,*)' 

£/«(-,!) 



(4.30) 



9n(±,*) 

we obtain the system of two recurrent relations 

Xn = fl{Xn-l,yn-l): (4.31) 

Vn = /2(x„_i,y„_i), 

with 

Pi{x,y) ^ a^ix^-' + b/i-iy'?-! + 2{Q - 1) 

^'^^'^^ R{x,y) fixi-^ + ^ji-^yi-^ + a + h + 2{Q-2y ^ ^ 

. P2{x,y) 6/2X^-1 + a/i-iy9-i + 2(Q - 1) 

f2{x,y) - 



R{x,y) + + a + 5 + 2(Q - 2) ' 

where the following notations are adopted 

a = exp{J + K), (4.33) 
b = exp(-J + Er), 
/i = exp(/i). 

In this approach statistical averages of all physical quantities can be expressed in terms of x and y variables 



defined by Eq. (4.31). For instance, when the total number of spins is n the magnetization along the first 



coordination axis which is the thermal average of the following form: 

= -E^^J'^) = -E^^^'^("-l))' (4.34) 

in terms of x and y are expressed as 



m 

n ^ — ' • - ■ fi 

i=l i=l 



ux'^ — u ^y'^ 

m = , r, (4.35) 

Hxi + fi-^yi + 2{Q-l)' ^ ' 

provided all sites of the lattice are equivalent. It is easy to see that this quantity plays the role of an order 
parameter of Ising universality class. Another order parameter belonging to the Q-state Potts universality 
class is 



1 " 

-Y.{6{a,A)), (4.36) 



p 

n 

which also can be regarded as a " quadrupolar" moment {S\ ). Thus 

= ^^'^ ^ ^"^^'^ f4 371 

^ ^^x1 + |JL-^y'i + 2{Q-ly ^ ' ' 

The order parameters m and "p define the three possible phases of the FCq model in case of the ferro- 
magnetic couplings, i.e. J > 0, ii' > 0, provided external magnetic field is vanished: 

(a) disordered (paramagnetic) phase: m = 0, p=l/Q 

(b) ferromagneticaly ordered phase: m 7^ , p 7^ l/Q 



(c) partially ordered (quadrupolar) phase: m = 0,p/l/Q 



4.4 Investigations of phase transitions in terms of dynamical systems 
theory 



In order to obtain physical results one must implement an iterative procedure for the RR ( [4.32 ). Namely, 
starting from the random initial conditions (xo,yo) one uses the simple iterations scheme and examines 



the behavior of physical quantities after a large number of iterations |156|. In a simplest case the iterative 
sequence {x„,y„} converges to a fixed point {x*,y*), which is defined by 



y* = f2{x*,y*) 



(4.38) 



This situation is inherent to the ferromagnetic case when both J and K are positive. As can easily be seen 
from Eqs. ( 4.34 ) - ( |4.37 ), when the magnetization and quadrupolar moment of the system take values m 
and p respectively, then the corresponding fixed point of the RR ( [4.32| ) is the following, provided h = 0: 



(Q-l)(p+m) 
1-p 

(Q-l)(p-m) 
1-p 



1/9 



1/g 



(4.39) 



Thus, according to the properties of different phases of the ferromagnetic FCq model one can establish the 
connection between them and the classification of the possible types of fixed points: 



(a) disordered (paramagnetic) phase: x* = y* = 1 

(b) ferromagneticaly ordered phase: x* / y* / 1 

(c) partially ordered (quadrupolar) phase: x* = y* ^ 1 



The condition x* = y* leads to (7„(+,l) = g„(— ,1), which means that the probability of spin "up" is 
equal to the probability of spin "down" which is really the case in paramagnetic phase. Applying the simple 
iterative scheme to Eqs. ( |4.31[ ), ( [4.35 ) and ( 4.37| ) one can obtain the plots of magnetization processes (m 
vs. H at fixed values of T) as well as the zero field magnetization and quadrupolar moment of the system. 

In Fig. 4.2 one can see the temperature dependencies of the order parameters of the FCg-model on 
Bethe lattice with coordination number q = 3 and Q = 2. As is obviously seen from Fig. [4.2| (a) in this case 
the zero field magnetization is always continuous, whereas the behavior of p is quite different. Depending 
on the value of ratio K/ J the p{T) curve could be continuous or can include one first order transition point. 
In Fig. |4.2| (b) one can see that at small values of K/ J only the one phase transition occurs from disordered 
phase to ferromagnetic phase omitting the partially ordered phase. This transition is of the second order 
until K/ J reaches some intermediate value above which one can see a discontinuity in the order parameter 
p at critical temperature corresponding to the transition between disordered and ferromagnetic phases. At 
large values of K/J the system undergoes successively two phase transitions at temperatures Tg{K) and 
Tf{K), the larger one corresponds to the transition between disordered and quandrupolar phases, the lower 
one - to the transitions from quandrupolar to ferromagnetic. Within this region of the values of K one can 
see the second order transition from disordered phase into "quadrupolar" phase and the further first order 
transition to the ferromagnetic phase. However, beginning with some value of K the transition between 
quadrupolar and ferromagnetic phases becomes continuous. 
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Figure 4.2: The temperature behavior of the order parameters for Q = 2 and q = 3. The transitions from disordered 
to ferromagnetic phase is of second order for all range of K .(a) Zero field magnetization for K/J = 0; 1 and 2.5; 

(b) Quadrupolar moment in zero field for K/J = 0; 1; 1.4 and 2. For K/J = 1 one can see the first order transition 
between disordered and ferromagnetic phases, whereas at K/J = 1.4 one can see two subsequent phase transitions, 
the continuous transition from disordered to quadrupolar phase and very close to that the first order transition to the 
ferromagnetic phase. With the further increase of K/ J the second transition becomes of the second order. 




Figure 4.3: The phase diagram of the model for Q — 2 and q — i. Solid line corresponds to the second order 
transitions, dotted - to the first order transitions. For small values of the ratio K/J one can see only single phase 
transition line between disordered hight-temperature phase and completely ordered ferromagnetic phase, whereas for 
K/J > 1.233 there is the partially ordered " quadrupolar" phase between them. The tricritical points are K/J w 
0.32; 1.233 and 1.6. 

In order to complete the picture of the phase structure of the model under consideration at Q = 2 we plot 
a phase diagram by separating the regions of the fixed points of different kinds of the RR from Eq. ( [4.32 ) in 



the {K/ T/ J)-plane. One can see in Fig. 4.3 that in case of Q = 2, q = 3 at X = the phase transition from 
disordered to ferromagnetic phase is of the second order. This feature maintains up to K/J ~ 0.32, where 
one can notice the tricritical point separating the region of continuous transitions from the region of the first 
order transitions. The line of the second order phase transitions between disordered and quadrupolar phases 
merges the line of transitions between disordered and ferromagnetic phase at K/J ~ 1.233, the latter, in its 
turn, again becomes of the second order at K/J « 1.6 and for large values of K/J it becomes parallel to 
the K/ J axis which means that the transition temperature between quadrupolar and ferromagnetic phases 
is unaffected by the value of K beginning with K/J ~ 4. 

The situation is quite different for Q > 2. In Fig. |4.4] the plots of zero field magnetization and quadrupo- 
lar moment for Q = 3 and q = 3 are presented. In this case both m and p are discontinuous at the critical 
temperatures. Moreover, this feature is preserved for p for all values of K/J, whereas m becomes continuous 
after K/J ~ 1.84. the corresponding phase diagram is presented in Fig. Here one can see a triple point 
at K/J ~ 1.67 where the merging of the three lines of the first order transitions takes place. There is also a 
tricritical point at K/J ~ 1.84 which separates the region of the first and the second order phase transitions 
on the line between quadrupolar and ferromagnetic phases. The line, as in case of Q = 2, goes to plateau 
beginning with K/J ~ 5, so the transition temperature again does not depend on K for K > 5 J. In general, 
for Q > 3 the pase diagrams of the model under consideration have the same topology as in case of Q = 3. 
The main feature is the decrease of the distance between triple and tricritical points with the increase of the 
Q. So, for Q = 30 one can obtain that this distance is of 10^^ order with K/J. Apparently, they have never 
shrunk to a single point at finite values of Q, but it could be the case when Q oo. The value of critical spin 





Figure 4.4: The temperature behavior of the order parameters for Q = 3 and g = 3. Here magnetization of the 
system undergoes a jump for small values of K/J and becomes continuous beginning with K/J w 1.84. Another order 
parameter always remains discontinuous at the transition points. In the plots p{T) for K/J =1.7 one can see that the 
transitions between " quadrupolar" and ferromagnetic phases is of the first order, whereas with the further increase of 
K/ J in becomes continuous. 



4,0 



3,5- 



3,0- 



2,5- 



I- 2,0- 



1,5- 



1,0- 



0,5- 




0,0 



0,5 



1,0 



1,5 



2,0 
K/J 



2,5 3,0 



3,5 



4,0 



Figure 4.5: The phase diagram of the model for Q = 3 and q = 3. Solid line corresponds to the second order 
transitions, dotted - to the first order transitions. The inset shows the merge point of the three first order lines 
(triple point): disordered-ferromagnetic, disordered-" quadrupolar" and "quadrupolar"-ferromagnetic. In this point 
at K/J pa 1.67 these three phase are in the equilibrium. The tricritical point on the line between "quadrupolar" and 
ferromagnetic phases is situated rather close to the triple point. The distance between them tends to zero with the 
increase of Q. 



component number Qc at which the phase transitions become of the first order at if = in the Face-cubic 
model on Bethe lattice is exactly equal to 3, which coincides with the mean-field results at oo, but, in 
contrast to that, on Bethe lattice this value is unaffected by the value of coordination number q. Obviously, 
this is the direct consequence of the fact, that the Bethe-Pierels approximation become exact on the Bethe 
lattice. 



Chapter 5 

Azimuthal Asymmetries in DIS as a 
Probe of Intrinsic Charm Content of the 
Proton i 



The notion of the intrinsic charm (IC) content of the proton has been introduced over 25 years ago in 
Refs [41, 42|. It was shown that, in the hght-cone Fock space picture |167, 168], it is natural to expect a 
five-quark state contribution to the proton wave function. The probabihty to find in a nucleon the five-quark 
component \uudcc) is of higher twist since it scales as 1/m^ where m is the c-quark mass [ |169(| . This com- 
ponent can be generated by gg cc fluctuations inside the proton where the gluons are coupled to different 
valence quarks. Since all of the quarks tend to travel coherently at same rapidity in the \uudcc) bound 
state, the heaviest constituents carry the largest momentum fraction. For this reason, one would expect 
that the intrinsic charm component to be dominate the c -quark production cross sections at sufficiently 
large Bjorken x. So, the original concept of the charm density in the proton [41, 42[ has nonperturbative 
nature and will be referred to in the present paper as nonperturbative IC. 

A decade ago another point of view on the charm content of the proton has been proposed in the 
framework of the variable flavor number scheme (VFNS) [^, The VFNS is an approach alternative to 
the traditional fixed fiavor number scheme (FFNS) where only light degrees of freedom (n, d, s and g) are 
considered as active. It is well known that a heavy quark production cross section contains potentially large 
logarithms of the type Ug In (Q^/m^) whose contribution dominates at high energies, — > oo. Within the 
VFNS, these mass logarithms are resummed through the all orders into a heavy quark density which evolves 



with according to the standard DGLAP [|170| , |171| , 172 [ evolution equation. Hence the VFN schemes 
introduce the parton distribution functions (PDFs) for the heavy quarks and change the number of active 
fiavors by one unit when a heavy quark threshold is crossed. We can say that the charm density arises 
within the VFNS perturbatively via the 5 — > cc evolution and will call it the perturbative IC. 

Presently, both perturbative and nonperturbative IC are widely used for a phenomenological descrip- 
tion of available data. (A recent review of the theory and experimental constraints on the charm quark 
distribution can be found in Refs. [173, 174]. See also Appendix 5.5 in the present thesis). In particular, 
practically all the recent versions of the CTEQ [^^ and MRST [^] sets of PDFs are based on the VFN 
schemes and contain a charm density. At the same time, the key question remains open: How to measure 
the intrinsic charm content of the proton? As mentioned in Introduction, the main theoretical problem is 
that production cross sections are not perturbatively stable. 



^The results considered in this chapter are published in Refs. 




Figure 5.1: Definition of tlie azimiitlial angle Lp in tlie micleon rest frame. 



In the present thesis we investigate the possibility to measure the charm content of the proton using the 
azimuthal asymmetries in heavy quark leptoproduction. For this purpose we calculate the IC contribution 
to azimuth-dependent process: 



+ N{p) ^ -q) + Q{pq) + X[Q]{px). 



(5.1) 



Neglecting the contribution of Z— boson as well as the target mass effects, the cross section of the reaction 
(13) unpolarized initial states may be written as 



1 



dxdQ2d(/5 (27r)2 xQ'^l-eV 



aT{x,Q^) + £aL{x,Q^) + e(7^(x, Q^) cos 2ip 
+2yje{l +£)ai{x, Q^) cos(^ 



(5.2) 



The quantity i measures the degree of the longitudinal polarization of the virtual photon in the Breit frame 

2(1 -y) 



and the kinematic variables are defined by 

S = 



+ P)\ 



y 



p_q 



i + {i-y?' 

= xyS, 



(5.3) 



X 



P 



2p-q' 

Am? 



(5.4) 



The cross sections cjj {i = T,L,A,I) in Eq. (|5.2| ) are related to the structure functions Fi{x,Q'^) as follows: 



Fi{x,Q^) 
F2{x,Q^) 



Sir'^aemX 



a,{x,Q'), {i = T,L,A,I) 



47r2 



Of, 



<T2(X,Q ), 



(5.5) 



where F2 = 2x{F'j' + F^) and (J2 = aj^ + a^. In Eq. (p^), (Tt {cfl) is the usual 7* A'" cross section describing 
heavy quark production by a transverse (longitudinal) virtual photon. The third cross section, a a, comes 
about from interference between transverse states and is responsible for the cos 2(/? asymmetry which occurs 
in real photoproduction using linearly polarized photons p^ , 35, 37 1. The fourth cross section, cj/, originates 



from interference between longitudinal and transverse components [ [L76| . In the nucleon rest frame, the 
azimuth ip is the angle between the lepton scattering plane and the heavy quark production plane, defined 
by the exchanged photon and the detected quark Q (see Fig. bA). The covariant definition of 99 is 



Qyi/x2+4myQ2 

cosLp = — - , sin 93= — = n • (5-6) 

y—r^y—nr 2y —r^y —nr 

= e'^-"V'?a^/3, n'^ = e^'^^^g.W/?. (5.7) 



In Eqs. (|5.4D and (|5.6|) , m and rufq are the masses of the heavy quark and the target, respectively. Usually, 
the azimuthal asymmetry associated with the cos 2^9 distribution, ^2i^(p> Q^)j is defined by 

^2(^(/0, a^, Q^) = 2(cos 2(/7)(p, Q^) 

d^cr;Ar((/9 = 0) + dVw((/? = 7r) - 2d^c7;Af(y? = 7r/2) 
dVw(v2 = 0) + dV/iv(y' = vr) + 2d3az7v(9' = ^/2) 

""^^"'"^'^ -A(x,Q^)" + ^f;'^'|, (5.8) 



f^T(x,Q2) +eaL(x,Q2) V ,^ ^l + ei?(x,Q2 

dxdQ2d(^ 



where d^cr;iv('/5) = (p, j;, Q2^ ^nd the mean value of cos rup is 



27r ^3^ 



/ d(/Jcosny? /\^'i^ {p,x,Q'^,Lp) 



{cos rup) {p, X, Q"^) = d^dQ dp ^ 

2^^ d^CTZTV 



if^^d^dg^^^'"'^''^) 



In Eq. (5^), the quantities R{x,Q ) and A{x,Q ) are defined as 



i?(x,Q2) = ^(^,Q2)^^(^^g2)^ (5^10) 

A(x,Q2) = ^(^^q2) ^2x^(x,Q2). (5.11) 

Likewise, we can define the azimuthal asymmetry associated with the cos p distribution, 
A^{p,x,Q^): 

A^{p,x,Q^) = 2{cosp) {p,x,Q^) (5.12) 

2dV;jv(y = 0) - 2d^aiN{p = vr) 

dVz7v(<^ = 0) + d^aiNip = vr) + 2dVi^(vp = 7r/2) 

2^6(1 + £) a/(x,Q^) _ 2^ / ^n , l + ^(x,Q^) 

" aT(x,g2)+eai(x,Q2) -^n^,^ '"^ ^ ''^ 1 + eR{x,Q^y 

where 

^/(x, Q2) = 2V2—{x, Q2) = 4^;r^(x, Q2). (5.13) 

Remember that y ^ 1 in most of the experimentally reachable kinematic range. Taking also into account 
that e = 1 + ©(y^), we find: 

A2^{p,x,Q^)=A{x,Q^) + 0{y^), A^{p,x,Q^) = Aj{x,Q^) + 0{y^). (5.14) 

So, like the ct2(x,(5^) cross section in the (/^-independent case, it is the parameters A{x,Q'^) and Aj{x,Q'^) 
that can effectively be measured in the azimuth-dependent production. 

In this paper we concentrate on the azimuthal asymmetry A{x, Q"^) associated with the cos 2(^-distribution. 
We have calculated the IC contribution to the asymmetry which is described at the parton level by the 



photon-quark scattering (QS) mechanism given in Fig. 5.2. Our main result can be formulated as follows: 



★ Contrary to the basic GF component, the IC mechanism is practically cos 2(/?- independent. This is 
due to the fact that the QS contribution to the £1^(2;, Q^) cross section is absent (for the kinematic 
reason) at LO and is negligibly small (of the order of 1%) at NLO. 



As to the 99-independent cross sections, our parton level calculations have been compared with the previous 
results for the IC contribution to a2{x, Q^) and 17^(2;, Q^) presented in Refs. |38, Apart from two trivial 



misprints uncovered in ||3^ for (Tl{x, Q^), a complete agreement between all the considered results is found. 

Since the GF and QS mechanisms have strongly different cos 299-distributions, we investigate the pos- 
sibility to discriminate between their contributions using the azimuthal asymmetry yl(2;,Q^). We analyze 
separately the nonperturbative IC in the framework of the FFNS and the perturbative IC within the VFNS. 

The following properties of the nonperturbative IC contribution to the azimuthal asymmetry within the 
FFNS are found: 

• The nonperturbative IC is practically invisible at low x, but affects essentially the GF predictions at 
large x. The dominance of the cos 2(/?-independent IC component at large x leads to a more rapid (in 
comparison with the GF predictions) decreasing of A(x, Q"^) with growth of x. 

• Contrary to the production cross sections, the cos 2(p asymmetry in charm azimuthal distributions is 
practically insensitive to radiative corrections at ~ m?. Perturbative stability of the combined 



GF+QS result for A{x, Q ) is mainly due to the cancellation of large NLO corrections in Eq. (5.11). 



• pQCD predictions for the cos 293 asymmetry are parametrically stable; the GF+QS contribution to 
A{x, Q^) is practically insensitive to most of the standard uncertainties in the QCD input parameters: 
HR, HF, J^QCD and PDFs. 

• Nonperturbative corrections to the charm azimuthal asymmetry due to the gluon transverse motion 
in the target are of the order of 20% at < m? and rapidly vanish at > m^. 

We conclude that the contributions of both GF and IC components to the cos2(^ asymmetry in charm 
leptoproduction are quantitatively well defined in the FFNS: they are stable, both parametrically and 
perturbatively, and insensitive (at > m?) to the gluon transverse motion in the proton. At large Bjorken 
X, the A{x^ Q^) asymmetry could be a sensitive probe of the nonperturbative IC. 



The perturbative IC has been considered within the VFNS proposed in Refs. |43, 44]. The following 
features of the azimuthal asymmetry should be emphasized: 

* Contrary to the nonperturbative IC component, the perturbative one is significant practically at all 
values of Bjorken x and > rm?. 

* The charm densities of the recent CTEQ and MRST sets of PDFs lead to a sizeable reduction (by 
about 1 /3) of the GF predictions for the cos 2(p asymmetry. 

We conclude that impact of the perturbative IC on the cos 2(p asymmetry is sizeable in the whole region of 
X and, for this reason, can easily be detected. 

Concerning the experimental aspects, azimuthal asymmetries in charm leptoproduction can, in principle, 
be measured in the COMPASS experiment at CERN, as well as in future studies at the proposed eRHIC 



1 43, m and LHeC ||9| colhders at BNL and CERN, correspondingly. 

The paper is organized as follows. In Section |5.1| we analyze the QS and GF parton level predictions 
for the (/3-dependent charm leptoproduction in the single-particle inclusive kinematics. In particular, we 
discuss our results for the NLO QS cross sections and compare them with available calculations. Hadron 
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Figure 5.2: The LO (a) and NLO (b and c) photon-quark scattering diagrams. 



level predictions for A{x,Q'^) are given in Section 5.2. We consider the IC contributions to the asymmetry 
within the FFNS and VFNS in a wide region of x and Q^. Some details of our calculations of the QS cross 



sections are presented in Appendix p.3[ An overview of the soft-gluon resummation for the photon-gluon 



fusion mechanism is given in Appendix 5.4. Some experimental facts in favor of the nonperturbative IC are 
briefly listed in Appendix |5.5|. 



5.1 Partonic Cross Sections 

5.1.1 Quark Scattering Mechanism 

The momentum assignment of the deep inelastic lepton-quark scattering will be denoted as 

m + Q{kQ) ^ m -q) + Q{pq) + X{px). (5.15) 
Taking into account the target mass effects, the corresponding partonic cross section can be written as 



follows [|4| 



dzdQ^V^ - (2vr)2 zQ- l-e L^^'«^"' " " '^^'^^^'^ 
+eaA,Q{z, A) cos 2(p + 2\/ e{l + e)ai^Q{z, A) cos 93 



(5.16) 



In Eq. (5.16|), we use the following definition of partonic kinematic variables: 



A 



m 



2 • 



In the massive case, the (virtual) photon polarization parameter, e, has the form |84| 

_ 2(l-y-AzV) 
^ l + (l-2/)2 + 2AzV' 



(5.17) 



(5.18) 



At leading order, 0{aem), the only quark scattering subprocess is 

7*iq) + QikQ)^QiPQ). 
The 7*(5 cross sections, o'^q {k = 2, L, A, I), corresponding to the Born diagram (see Fig. 

aB{z)VlT4X^ 6{1 - z), 



(5.19) 



a) are: 



(0) 



a 



(0) 
A,Q 



(z,A) 



0, 



6(1 



(5.20) 



with 



0-5(2:) 



a(°)(.,A) 



(5.21) 



where cq is the quark charge in units of electromagnetic coupling constant. 

To take into account the NLO ©(aemOs) contributions, one needs to calculate the virtual corrections to 
the Born process (given in Fig. |5.2| c) as well as the real gluon emission (see Fig. |5.2| b): 

l*{q) + Q{kQ)^Q{pQ)+9{Pg). (5.22) 

The NLO (/9-dependent cross sections, (T^q and (t^^q, are described by the real gluon emission only. 
Corresponding contributions are free of any type of singularities and the quantities <t^q and irj^Q can be 
calculated directly in four dimensions. 

In the (^-independent case, a^2Q ^^'^ "^lq' work in four dimensions. The virtual contribu- 



tion (Fig. |5.2| c) contains ultraviolet (UV) singularity that is removed using the on-mass-shell regularization 
scheme. In particular, we calculate the absorptive part of the Feynman diagram which has no UV diver- 
gences. The real part is then obtained by using the appropriate dispersion relations. As to the infrared (IR) 
singularity, it is regularized with the help of an infinitesimal gluon mass. This IR divergence is cancelled 
when we add the bremsstrahlung contribution (Fig. ^.2f 3). Some details of our calculations are given in 
Appendix |5.3| . 

The final (real+virtual) results for j*Q cross sections can be cast into the following form: 

^C7^aB(l)Vl + 4A(5(l -z)|-2 + 41nA- V1 + 4A Inr 
1 + 2A r 



+ 



V1 + 4A 



2Li2(r2) +4Li2(-r) + 31n^(r) - 4 Inr 
+4 In r ln(l + 4A) - 2 In r In Aj | 



+ 4^^^^^^"\l+4Az2)3/2|[i_(i_;,),]2 



1 



1 - 3z - iz^ 



+6z^ + 8z^ - 8z' 



+6Xz (3 - 18z + 13z^ + lOz^ 



+4X^z'^ (8 - 77z + 65z^ - 2/) 



+16X^z^ (1 - 21z + 12^2) - l28A^z^ 
(1 + 2 + 2^2 + 2z^) + 2Az (2 - llz - 11^2) 

+8X'^z'^ (1 - 9z) 

2^4 



(5.23) 



+ 



21n7:'(z,A) 
Vl + 4Az2 L 



1 + 4Xyz^ 4(1 + 2A)(1 + 4A)2z4 lnD{z, A) 



(1-^)- 



Vl + 4A2 



(1 



>,A) = ^CFaB{l)^j^t=5{l-z)\-2 + 4lnX 



+ 



4A 



Inr + 



1 + 2A 



, .... , , 2Li2(r2) +4Li2(-r) 

+31n2(r) - 41nr + 41nr ln(l + 4A) - 21nrlnA | 

-\z (13 - 19z - 2z^ + 8z^) 
-2\^z^ (31 - 39z + 8^2) 
-SA^z^ (10 - 7z) - 32A^z^ 

[3 + 3z + 16Az] 



2Az^ \nD{z,\) 
Vl + 4Az2 

8A(1 + 4A)z4 4A(1 + 2A)(1 + 4.\)z^ lnZ?(z, A) 



(1 - 



Vl + 4Az2 (l-^)_ 



(5.24) 



z{l-z) 



2tx 



(l + 4Az2)3/2 (i_ A)z 

2AlnL>(z,A) 
7l + 4Az2 



[1 + 2A(4 - 3z) + 8A2 



+- 



[2 + z + 4Az] L 



(5.25) 



^i;^(z,A) 



8^/2 



Cf&b{z) 



1 



(l+4Az2)2 [1 _ (1 _ ;^)^]3/2 



-(l-z)(l + 2z) 



4Az (10 - lOz - z^ + 2z^ 



-8A^z2 (25 - 29z + 82^) - 96A^z^ (3 - 2z) - 128A 
+8VA2 [1 - (1 - A)z] [1 - z^ + Az(13 - llz) + 4A2z2(7 - 4z) + leA^z^] |. 
In Eqs. (|5.23H5.26| ), = (N^ - l)/(2iVc), where iVc is number of colors, while 



B{z,X) 



1 + 2Az - Vl + 4Az2 



l + 2Az + Vl + 4Az2' 
The so-called "plus" distributions are defined by 



r=^D{z = \,\) 



yi + 4A - 1 
VI + 4A + 1 ' 



(5.26) 



(5.27) 



\g{z)\^=g{z)-b(^-z) J dCgiC). 





For any sufficiently regular test function h{z), Eq. ( 5.2^ ) gives 



1-Z 



, In^(l-z) „ , , , , , , ,ln*^+i(l-a) 
dz^ ^-[h{z) -h{l)]+h{l)- ^ ' 



l-z 



k + l 



(5.28) 



(5.29) 



To perform a numerical investigation of the inclusive partonic cross sections, (T^^q (A; = T, L, A, I), it is 



convenient to introduce the dimensionless coefficient functions c 

el)aemasiP''^'^ ~ 



(n,0 
k,Q ' 

n 



0"fc,Q(??,A,/X ) 



5:(4vra.(^2))n^^(n.)(^^A)l,, 



n=0 



/=0 



(5.30) 




where /x is a factorization scale (we use ^ = = f^R) and the variable r] measures the distance to the 
partonic threshold: 

r?=4-l = lf^, s={q + kQ)\ (5.31) 



Our analysis of the quantity c^'q (^/, A) is given in Fig. 5.3. One can see that c^'g is negative at low 
(A^^ < 1) and positive at high (A^^ > 20). For the intermediate values of Q^, c^'g (?7, A) is an alternating 
function of rj. 

.(0,0), 



Our results for the coefficient function Cjq {r],X) at several values of A are presented in Fig. 5^. It 
jeen tl 
function: 



is seen that cf''Q is negative at all values of r] and A. Note also the threshold behavior of the coefficient 



c 



m^rj - 0, A) = -Vln'Cp-^ + 0(r?). (5.32) 



This quantity takes its minimum value at Am = 1/4: 

. Let us analyze the numerical significance of the cosip- and cos 299-distributions for the QS component. 
It is difficult to compare directly the (T^g(z,A) and a^^Q{z,X) cross sections given by the usual functions 
( ^.25 ) and ( ^.261) with the (/3-independent contributions ^^q{z, A) and ^^q{z, A) described by the generalized 
functions ( |5.2[1| ) and ( 5.23] ). For this reason, we consider the Mellin moments of the corresponding quantities 
defined as 



ai,Q{N,X) = J ai^Q{z,X)z''-'dz, {i = 2,L,A,I). (5.33) 



The Mellin transform of the Born level cross sections is trivial: a (°^(Ar,A) = aB(l)Vl + 4A. The Mellin 
moments of the NLO results have been calculated numerically. We use for OsifJ-p) the one- loop approximation 
with A4 = 326 MeV, hf = + Q2 and m = 1.3 GeV. 

The left panel of Fig. U presents the ratio o^a,q^^ ^ A)/4°^(iV, A) as a function of for several values 
of variable A: A'^ = 1,4,10,20 and 100. One can see that this ratio is negligibly small (of the order of 1%). 
Moreover, our analysis shows that the ratio cJ^g(A^, A)/(T2'g(A^, A) is less than 1.5% for all values of A and 

> 0. This implies that the photon-quark scattering contribution is practically cos 2(/?-independent and, 
for this reason, we will neglect (7a,q{_z^X) cross section in our further analysis. 

In the right panel of Fig. 5^, the A^-dependence of the ratio 2\plaf^Q{N , X)/a^Q{N, A) is given for the 
same values of A. One can see that this ratio is of the order of 10-15% at small N and sufficiently high Q^. 
This fact indicates that the cos 99-distribution caused by the QS component may be sizable. 
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Figure 5.4: The quantities o-^^.q(^, A)/CT^°^(iV, A) (left panel) and 2^/2 ir}^^ (iV, A)/(T^°^(7V, A) {right panel) at several 
values of A. 



5.1.2 Comparison with Available Results 

For the first time, the NLO 0{a 

emCts) corrections to the 93-independent IC contribution have been calculated 
a long time ago by Hoffmann and Moore (HM) [^]. However, authors of Ref. [^8| don't give explicitly their 
definition of the partonic cross sections that leads to a confusion in interpretation of the original HM results. 
To clarify the situation, we need first to derive the relation between the lepton-quark DIS cross section, daig, 
and the partonic cross sections, fj^^^ and a^^\ used in [^]. Using Eqs. (C.l) and (C.5) in Ref. [^], one 
can express the HM tensor aj^ in terms of "our" cross sections a2,Q and o'l^q defined by Eq. ( 5.16| ) in the 
present paper. Comparing the obtained results with the corresponding definition of aj^ via the HM cross 
sections a^^^ and a^^^ (given by Eqs. (C.16) and (C.17) in Ref. we find that 



^l,q{z,X) 



aB{z)VT+4Xz^a^^\z,X), 



2a Biz) 



(5.34) 
(5.35) 



VI + 4Az2 

Now we are able to compare our results with original HM ones. It is easy to see that the LO cross sections 



(defined by Eqs. (37) in |Q and Eqs. ( ^.201) in our paper) obey both above identities. Comparing with 
each other the quantities cij^^ and a^Q (given by Eq. (51) in [^] and Eq. ( ^.23 ) in this paper, respectively). 



we find that identity (5.34) is satisfied at NLO too. The situation with longitudinal cross sections is more 
complicated. We have uncovered two misprints in the NLO expression for a^^^ given by Eq. (52) in p^]. 
First, the r.h.s. of this Eq. must be multiplied by z. Second, the sign in front of the last term (proportional 
to (5(1 — z)) in Eq. (52) in Ref. [38| must be changed 0. Taking into account these typos, we find that relation 
( ^.35 ) holds at NLO as well. So, our calculations of (T2,q and o'l,q agree with the HM results. 

Recently, the heavy quark initiated contributions to the independent DIS structure functions, F2 and 
Fl, have been calculated by Kretzer and Schienbein (KS) (3^. The final KS results are expressed in terms 
of the parton level structure functions Hf and /^2- Using the definition of Hf and H2 given by Eqs. (7, 8) 
in Ref. [l39|, we obtain that 



'72,q(^,A) 



^T,Qiz,X) 

2tt 



a. 



asiz) Hf{e,X) 



2tt VI + 4A VI + 4Az2 ' 

aB{z)^ 



(5.36) 



1 + 4A 
1 + 4Az2 



^l(e',A), 



■^Note that this term originates from virtual corrections and the virtual part of the longitudinal cross section given by Eq. (39) 
in Ref. pa^ also has wrong sign. See Appendix 5.3 for more details. 




Figure 5.5: The LO photon-gluon fusion diagrams. 



where (Jt,q = ^2,Q — ^l,q and ^l,q are defined by Eq. ( 5.16 ) in our paper and ^' = z (l + vT-i-^) / 
(l + Vl + 4Az2^ . To test identities (|]36|), one needs only to rewrite the NLO expressions for the functions 
Hf{£^', A) and -fff ('^'' (given in Appendix C in Ref. [^) in terms of variables z and A. Our analysis shows 
that relations (|5.36|) hold at both LO and NLO. Hence we coincide with the KS predictions for the 7*Q 
cross sections. 



However, we disagree with the conclusion of Ref. |39| that there are errors in the NLO expression for 
cj(2) given in Ref. |38| ^. As explained above, a correct interpretation of the quantities a^"^^ and a^^^ used in 
leads to a complete agreement between the HM, KS and our results for 99-independent cross sections. 

As to the (/^-dependent DIS, pQCD predictions for the 7*Q cross sections aA,Q{z, A) and ai^Q{z, A) in the 
case of arbitrary values of m? and are not, to our knowledge, available in the literature. For this reason, 
we have performed several cross checks of our results against well known calculations in two limits: m? 
and —>■ 0. In particular, in the chiral limit, we reproduce the original results of Georgi and Politzer 
|177| and Mendez |178| ] for aj^Q{z, A — > 0) and aA,Q{z, A ^ 0). In the case of —>■ 0, our predictions for 
^2,q{s,Q'^ — > 0) and aA,Q{s,Q'^ — > 0) given by Eqs. ( p.23| , 5?25| ) reduce to the QED textbook results for the 
Compton scattering of polarized photons ]17S|]. 



5.1.3 Photon-Gluon Fusion 

The gluon fusion component of the semi-inclusive DIS is the following parton level interaction: 

l{e) + g{kg) ^ l{i -q) + Q{pq) + X[Q]{px). 



(5.37) 



Corresponding lepton-gluon cross section, d&ig, has the following decomposition in terms of the helicity 'y*g 
cross sections: 



Or 



1 



dzdQMip (27r)2 zQ^ 1 - e 



^2,g{z, A) - (1 - e)aL,g{z, A) 



+eaA,g{z, A) cos 2ip + e(l + e)ai^g{z, A) cos (p 



(5.38) 



where the quantity e is defined by Eq. ( [5.3[ ) with y = {q ■ kg)/ {i ■ kg). 
At LO, 0{ 

Q^emQ^s)) the only gluon fusion subprocess responsible for heavy flavor production is 



l*{q) + 9{kg) ^ Q{pq) + Qip^). 



(5.39) 



The 7*5 cross sections, {k = 2, L, A, /), corresponding to the Born diagrams given in Fig. |5.5| have the 



^In detail, the KS point of view on the HM resuhs is presented in PhD thesis |^^, pp. 158-160. 



form pO, pi: 



(0) 



2-K 



[(1 - zf + + 4Az(l - 3z) - SX^z^] In - 



1 + 



vr 



(z)z{-2Azlni±|i + (l 



[l + 4z(l-z)(A-2)]/3,}, 



A,9 



a 



(0) 



;^,A) 



— ctb(z)z<^ 2A [1 - 2z(l + A)] In - 
vr L 1 



1 + 



/3. 



+ (1 



2A)(l-z)/3,}, 



0, 



where ^^(z) is defined by Eq. ( |5.21 ) and the following notations are used: 



2q ■ kg 



m 



I3z 



4Xz 



(5.40) 

(5.41) 

(5.42) 
(5.43) 



Note that the cos if dependence vanishes in the GF mechanism due to the Q 



(5.44) 

Q symmetry which, at 
leading order, requires invariance under ip ^ ip + tt [182]. 

As to the NLO results, presently, only (/^-independent quantities o-g"'^^ and o"^^^ are known exactly [p^ ]. 
For this reason, we will use in our analysis the so-called soft-gluon approximation for the NLO 'y*g cross 



sections (see Appendix |5.4|) . As shown in Refs. |3^, |3^, at energies not so far from the production 
threshold, the soft-gluon radiation is the dominant perturbative mechanism in the 'y*g interactions. 



5.2 Hadron Level Results 

5.2.1 Fixed Flavor Number Scheme and Nonperturbative Intrinsic Charm 

In the fixed flavor number scheme the wave function of the proton consists of light quarks u, d, s and 
gluons g. Heavy flavor production in DIS is dominated by the gluon fusion mechanism. Corresponding 
hadron level cross sections, ^"^^^'^(x, A), have the form 

1 

crk,GFix,X) = J dzg{z,fiF)^k,gix/z,X,iJ,F) , {k = 2,L,A,I), (5.45) 

X 

X = x(l + 4A), (5.46) 

where g{z,fiF) describes gluon density in the proton evaluated at a factorization scale jip. The lowest 
order GF cross sections, a^^^ {k = 2,L,A,I), are given by Eqs. ( ^.411 ). The NLO results, to the 

next-to-leading logarithmic accuracy are presented in Appendix |5. 4 

We neglect the ^*q{q) fusion subprocesses. This is justified as their contributions to heavy quark lepto- 



production vanish at LO and are small at NLO |23]. 



In the FENS, the intrinsic heavy flavor component of the proton wave function is generated by gg QQ 
fluctuations where the gluons are coupled to different valence quarks. In the present paper, this component 
is referred to as the nonperturbative intrinsic charm (bottom). The probability of the corresponding five- 
quark Fock state, \uudQQ^, is of higher twist since it scales as Ag^^^/m^ |16£]. However, since all of the 



quarks tend to travel coherently at same rapidity in the \uudQQ) bound state, the heaviest constituents 
carry the largest momentum fraction. For this reason, the heavy flavor distribution function has a more 



*This approach is sometimes referred to as the fixed-order perturbation theory (FOPT). 



"hard" z-behavior than the Ught parton densities. Since ah of the densities vanish at z — > 1, the hardest 
PDF becomes dominant at sufficiently large z independently of normalization. 

Convolution of PDFs with partonic cross sections does not violate this observation. In particular, 
assuming a gluon density g(^z) ~ (1 — z)" (where n = 3 — 5), we obtain that the LO GF contribution to 
scales as (1 — x)""*"^^^ at x ^ 1) where x is defined by Eq. (5.46). In the case of Hoffman and Moore charm 



density (see below), the LO IC contribution is proportional to (1 — x) at x ~^ 1- It is easy to see that, 
independently of normalizations, the IC contribution to be dominate over the more "soft" GF component 
at large enough x. 

For the first time, the intrinsic charm momentum distribution in the five-quark state \uudcc) was derived 
by Brodsky, Hoyer, Peterson and Sakai (BHPS) in the framework of a light-cone model |4^, Neglecting 
the transverse motion of constituents, they have obtained in the heavy quark limit that 

c{z) = —z^ \6z{l + z)lnz + {l-z){l + Wz + z'^)] , (5.47) 
6 

where = 36 corresponds to a 1% probability for IC in the nucleon: c{z)dz = 0.01. 

Hoffmann and Moore (HM) [^8| incorporated mass effects in the BHPS approach. They first introduced 
a mass scaling variable ^, 

e= ^^^=, a = i±^^, (5.48) 

where Xn = /Q'^ ■ To provide correct threshold behavior of the charm density, the constraint ^ < 7 < 1 
was imposed where 

7 = , x = - — . (5.49) 

1 + 1 + 4A-A^ 

Resulting charm distribution function, 0(^,7), has the following form in the HM approach: 

c{(,l) = { 7 (5.50) 

0, e>7 



with c(^) defined by Eq. (5.47). Corresponding hadron level cross sections for the (c + c) production, 



Ck^Qsix, X), due to the heavy quark scattering (QS) mechanism, are 



7 

f dz 

crk,Qsix,X)= / c+{z,-/)ak,c{C/z,X) , {k = 2,L,A,I), (5.51) 

J y' 1 + 4Xi,^ / z'^ 



where the charm density c+(z,7) = 0(2,7) + 0(2;, 7). The LO and NLO expressions for the partonic cross 
sections (Ji:^c{z,X) are given by Eqs. ( p. 201 ) and ( |5.23| - |5.2^ ), respectively. 

Note also that, in the FENS, the full cross section for the charm production, ak{x,X), is simply a sum 
of the GF and IC components: 

ak{x,X) = ak,GFix,X) + ak,Qsix,X), {k = 2,L,A,I). (5.52) 

Let us discuss the FENS predictions for the hadron level asymmetry parameter A{x,Q'^) defined by 
Eq. ( |5TlD . Fig. U shows A{x, X) function of x for several values of variable A: A ^ = 1,4, 10, 20, 50 



and 100. We display both LO and NLO predictions of the GF mechanism as well as the analogous results of 
the combined GF+QS contribution. The azimuthal asymmetry due to the mere LO GF component is given 
by solid line. The NLO GF predictions are plotted by dashed line. The LO and NLO results of the total 
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Figure 5.6: Azimuthal asymmetry parameter A{x, A) in tlie FFNS at several values of A in the case 

The following contributions are plotted: GF^^^^ (solid lines), GF^^^'+fcr-kick (dotted lines), GF^'^'^^^ (dashed lines), 

GF(lo)_^qs(lo) (dash-dotted lines) and GF^'^^^^+QS^'^^^) (long-dashed lines). 



GF+QS contribution are given by dash-dotted and long-dashed hnes, respectively. In our calculations, the 



CTEQ5M [183] parametrization of the gluon distribution function is used and a 1% probability for IC in 
the nucleon is assumed. Throughout this paper, the value jip = fJ-R = \/ fn^ + for both factorization and 
renormalization scales is chosen. In accordance with the CTEQ5M parametrization, we use rric = 1.3 GeV 
and A4 = 326 MeV p^ ]. 

One can see from Fig. 5.6 the following basic features of the azimuthal asymmetry, A{x,\), within the 
FFNS. First, as expected, the nonperturbative IC contribution is practically invisible at low x, but affects 
essentially the GF predictions at large x. Since, contrary to the GF mechanism, the QS component is 
practically cos 2(^- independent, the dominance of the IC contribution at large x leads to a more rapid (in 
comparison with the GF predictions) decreasing of A) with growth of x. 

The most remarkable property of the azimuthal asymmetry is its perturbative stability. In Refs. 
the NLO soft-gluon corrections to the GF predictions for the cos 293 asymmetry in heavy quark photo- and 
leptoproduction was calculated. It was shown that, contrary to the production cross sections, the quantity 
A{x, A) is practically insensitive to soft radiation. One can see from Fig. in the present paper that 
the NLO corrections to the LO GF predictions for A[x, A) are about few percent at not large x. This 
implies that large soft-gluon corrections to ct^q]? and cr^Qp (increasing both cross sections by a factor 
of two) cancel each other in the ratio (<^^f(;F (•^' ■^) '^ith a good accuracy. In terms of so-called 

if-factors, Kk{x,X) = (u^.^^'^ycr^^'^^) (x. A) for k = 2,L,A,I, perturbative stability of the GF predictions 
for A{x, A) is provided by the fact that the corresponding /C-factors are approximately the same at not large 
x: Ka,gf{x,X) ^ K2,gf{x,X)- 

Comparing with each other the dash-dotted and long-dashed curves in Fig. 5.6, we see that the NLO 
corrections to the combined GF+QS result for A{x, A) are also small. In this case, three reasons are 
responsible for the closeness of the LO and NLO predictions. At small x, where the nonperturbative IC 
contribution is negligible, perturbative stability of the asymmetry is provided by the GF component. In 
the large-x region, where the IC mechanism dominates, the azimuthal asymmetry rapidly vanishes with 
growth of X at both LO and NLO because the QS component is practically cos 2(/7- independent, cr^^S.^'' ■^) ~ 
^^Ac^'^^ = Q- intermediate values of x, where both mechanisms are essential, perturbative stability 
of A{x, A) is due to the similarity of the corresponding K-factors: K2,gf{x, A) ~ K2,qs{x, A) ^. 

Another remarkable property of the azimuthal asymmetry closely related to fast perturbative convergence 
is its parametric stability Q The analysis of Refs. |34, ^ shows that the GF predictions for the cos2(^ 
asymmetry are less sensitive to standard uncertainties in the QCD input parameters (m, fiR, fip, -^qcd and 
PDFs) than the corresponding ones for the production cross sections. We have verified that the same 
situation takes also place for the combined GF+QS results. 

Let us discuss how the GF predictions for the azimuthal asymmetry are affected by nonperturbative 
contributions due to the intrinsic transverse motion of the gluon in the target. Because of the relatively 
low c-quark mass, these contributions are especially important in the description of the cross sections for 
charmed particle production. 

To introduce degrees of freedom, kg ~ + ^T, one extends the integral over the parton distribution 



^Although the ratio (A'^^'^^/^*-^'^^)(a;, A) is sizeable at sufficiently large x, the absolute values of the quantities j4'^'^'(x, A) 
and A'^^"^' (a;, A) become so small that it seems reasonable to consider the asymmetry as equally negligible at both LO and 
NLO and treat the predictions as perturbatively stable. 

^Note however that this similarity takes only place at intermediate values of x where both GF and QS components are 
essential. In the low- and large- a; regions, the factors K2,gf{x, X) and K2,qs{x, X) are strongly different. 

^Of course, parametric stability of the fixed order results does not imply a fast convergence of the corresponding series. 
However, a fast convergent series must be parametrically stable. In particular, it must be hr- and /iF-independent. 



function in Eq. ( |5.45| ) to kx-space, 

dCgiCf^p) ^ dC d^kTf{kT)g{C,^^F)■ (5.53) 

The transverse momentum distribution, fikx), is usuahy taken to be a Gaussian: 



In practice, an analytic treatment of fcy effects is usually used. According to Ref. |184], the /cy-smeared 
differential cross section of the process (^]^) is a two-dimensional convolution: 

dxdQ^PQTd^ (pV) = /d'^T^^^^^^^,^^(pV - \kT)- (5.55) 



The factor ^ in front of fey in the r.h.s. of Eq. ( 5. 551) reflects the fact that the heavy quark carries away 



about one half of the initial energy in the reaction ( |5.1| ). 

Values of the fc^-kick corrections to the LO GF predictions for the cos 2ip asymmetry in the charm 
production are shown in Fig. 5.(; by dotted curves. Calculating the /c-r-kick effect we use {k"^) = 0.5 GeV'^. 
One can see that fcr-smearing for ^(x,Q2) is about 20-25% in the region of low < vn? and rapidly 
decreases at high Q^. 



In Fig. [5.7| , the dependence of the asymmetry A{x, A) on the nonperturbative intrinsic charm content of 
the proton is presented. We plot the LO predictions for A{x^ A) as a function of x for several values of the 
variable A and quantity Pc = c{z)dz describing a probability for IC in the nucleon. Dash-dotted curves 
describe the GF^^^^-^QS^^^) contributions with = 5%, 1%, 0.1% and 0.01%. Solid lines correspond 



to the case when Pc = 0. Comparing with each other Figs. 5.6 and 5.7, one can see that even a 0.1% 
contribution of the nonperturbative IC to the proton wave function could be extracted from the cos 2ip 
asymmetry at large enough Bjorken x. 

5.2.2 Variable Flavor Number Scheme and Perturbative Intrinsic Charm 



One can see from Eqs. ( ^.41 ) that the GF cross section (72^(2, A) contains potentially large logarithm, 



ln((5^/m^). The same situation takes also place for the QS cross section al2^Q{z,\) given by Eq. (5.23). At 
high energies, — > 00, the terms of the form asln((5^/m^) dominate the production cross sections. To 
improve the convergence of the perturbative series at high energies, the so-called variable flavor number 
schemes (VFNS) have been proposed. Originally, this approach was formulated by Aivazis, Collins, Olness 



and Tung (ACOT) |18|, ^|. 



In the VFNS, the mass logarithms of the type q" ln"'((5^/m^) are resummed via the renormalization 
group equations. In practice, the resummation procedure consists of two steps. First, the mass logarithms 
have to be subtracted from the fixed order predictions for the partonic cross sections in such a way that in 
the asymptotic limit ^ 00 the well known massless MS coefficient functions are recovered. Instead, a 
charm parton density in the hadron, c(x,(5^), has to be introduced. This density obeys the usual massless 
NLO DGLAP evolution equation with the boundary condition c{x,Q^ = Qq) = where Qq ~ m^. So, we 
may say that, within the VFNS, the charm density arises perturbatively from the (7 — > cc evolution. 

In the VFNS, the treatment of the charm depends on the values chosen for Q^. At low < Qq, 
the production cross sections are described by the light parton contributions {u,d,s and g). The charm 
production is dominated by the GF process and its higher order QCD corrections. At high ^ m?, the 
charm is treated in the same way as the other light quarks and it is represented by a charm parton density 




Figure 5.7: The LO predictions for A{x, A) in the FFNS at several values of A and Pc = /q c(z)dz. Dash-dotted 
curves describe the GF^^^'+QS^^^^ contributions with = 5%, 1%, 0.1% and 0.01%. Solid lines correspond to the 
case when Pc = 0. 



in the hadron, which evolves in . In the intermediate scale region, ~ m^, one has to make a smooth 
connection between the two different prescriptions. 

Strictly speaking, the perturbative charm density is well defined at high ^ but does not have a 
clean interpretation at low . Since the perturbative IC originates from resummation of the mass logarithms 
of the type a" ln"'(Q^/m^), it is usually assumed that the corresponding PDF vanishes with these logarithms, 
i.e. for Q'^ < Q\k, tit?. On the other hand, the threshold constraint W'^ = {q + p)"^ = Q'^{l/x — 1) > 4m^ 
implies that Qq is not a constant but "live" function of x. To avoid this problem, several solutions have been 



proposed (see e.g. Refs. |186, |188| ). In this paper, we use the so-called ACOT(x) prescription [ |186| | which 
guarantees (at least at > m?) the correct threshold behavior of the heavy-quark-initiated contributions. 
Within the VFNS, the 99- independent charm production cross sections have three pieces: 

cj2(x, A) = a2,GFix, A) - a2,suB{x, A) + cr2,Qs{x, A), (5.56) 

where the first and third terms on the right hand side describe the usual (unsubtracted) GF and QS 
contributions while the second (subtraction) term renders the total result infra-red safe in the limit — > 0. 
The only constraint imposed on the subtraction term is to reproduce at high energies the familiar MS 
partonic cross section: 

lim [<T2,g(z, A) - (T2,sub{z, A)] = affiz). (5.57) 

Evidently, there is some freedom in the choice of finite mass terms of the form A" (with a positive n) in 
(T2,5(7b(-z, A). For this reason, several prescriptions have been proposed to fix the subtraction term. As 



mentioned above, we use the so-called ACOT(x) scheme (|186|| . 



According to the ACOT(x) prescription, the lowest order (/3-independent cross section is 



1 

2 



a'ilix/z, A) - ^ In ^ as {x/z) {x/z) 



(5.58) 



X 



+aBix)c+{x,fJ'F) 



where Pg'^c is the LO gluon-quark sphtting function, Pg-idC) = [(1 — C)^ + C^] / 2, and the LO GF cross 



section is given by Eqs. ( 5.41 ). Remember also that x = 2^(1 + 4A) and c+(C,/^f) = c{C, ^f) +c(C,/Uf)- 



The asymptotic behavior of the subtraction terms is fixed by the parton level factorization theorem. This 
theorem implies that the partonic cross sections da can be factorized into process-dependent infra-red safe 
hard scattering cross sections da, which are finite in the limit m — > 0, and universal (process-independent) 
partonic PDFs fa^i and fragmentation functions dn^Q- 

da{j* + a^Q + X)=Y, fa-.i{0 ^ da(7* + i ^ n + X) <g) dn^qiz). (5.59) 

i.n 



In Eq. ( 5.591 ), the symbol denotes the usual convolution integral, the indices a, z,n and Q denote partons. 



Pi = Qpa and pq = zpn- All the logarithms of the heavy-quark mass (i.e., the singularities in the limit 
m — > 0) are contained in the PDFs fa-*i and fragmentation functions dn~*Q while da are IR-safe (i.e., are 
free of the In terms) . The expansion of Eq. ( [5. 591) can be used to determine order by order the subtraction 



terms. In particular, for the LO GF contribution to the charm leptoproduction one finds |43] 

[z, In it^l/m^)) = f^^l, (C, In (mfM')) ® ^khs(^/0, (k = 2, L, A, /), (5.60) 



^k.SUB 



where /^Hc (C; In (/i|,/m^)) = (as/27r) In (^|^/m^) Pg!?lc (C) describes the charm distribution in the gluon 
within the MS factorization scheme. 
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Figure 5.8: Azimuthal asymmetry parameter A{x, A) in tlie VFNS at several values of A. Tlie following contributions 
are plotted: GF^^°^ (solid curves), GF^^^^-SUB^^^^+QS^^^^ with the CTEQ5M set of PDFs (dotted curves) and 
Qp(LO)_guB^^°^ QS(LO) ^i^jj ^Y^^ CTEQ4M set of PDFs (dashed curves). 



One can see from Eq. ( 5.60| ) that the azimuth-dependent GF cross sections aA,GF and ai^cF don't have 



subtraction terms at LO because the lowest order QS contribution is 99- independent. For this reason, the 
cos 299-dependence within the VFNS has the same form as in the FFNS: 

1 

a^j^'^\x,\)= / dzg{z,fiF)cF'A,„{x/z,X) . (5.61) 



Fig. |5.8| shows the ACOT(x) predictions for the asymmetry parameter A{x, A) at several values of variable 
A: A^^ = 1, 4, 10, 20, 50 and 100. For comparison, we plot also the LO GF predictions (solid curves). In the 
ACOT(x) case, we consider the CTEQ5M (dotted lines) and CTEQ4M (dashed curves) parametrizations of 
the gluon and charm densities in the proton. Corresponding values of the charm quark mass are = 1.3 



GeV [HI (for the CTEQ5M PDFs) and = 1.6 GeV (for the CTEQ4M PDFs). The default value 
of the factorization scale is hf = -s/ m? + Q^. 



One can see from Fig. the following properties of the azimuthal asymmetry, A{x,X), within the 
VFNS. Contrary to the nonperturbative IC component, the perturbative one is significant practically at all 
values of Bjorken x and > m?. The perturbative charm contribution leads to a sizeable decreasing of 
the GF predictions for the cos 2(^-asymmetry. In the ACOT(x) scheme, the IC contribution reduces the 
GF results for j4(x. A) by about 30%. The origin of this reduction is straightforward: the QS component is 
practically cos 2</?-independent. 

The ACOT(x) predictions for the asymmetry depend weakly on the parton distribution functions we 
use. It is seen from Fig. |5.8| that the CTEQ5M and CTEQ4M sets of PDFs lead to very similar results for 
A{x^ A). Note that we give the CTEQ5M predictions at low x only because of irregularities in the CTEQ5M 
charm density at large x. 

We have also analyzed how the VFNS predictions depend on the choice of subtraction prescription. 



In particular, the schemes proposed in Refs. [39, |188U have been considered. We find that, sufficiently 
above the production threshold, these subtraction prescriptions reduce the GF results for the asymmetry 
by approximately 30 50%. 

One can conclude that impact of the perturbative IC on the cos 2(p asymmetry is essential in the whole 
region of Bjorken x and therefore can be tested experimentally. 

5.3 Appendix A: Virtual and Soft Contributions to the Quark Scattering 

In this Appendix we reproduce some results of Hoffmann and Moore for the (/j-independent QS cross sections, 
and correct two misprints uncovered in Ref. [^]. We work in four dimensions, in the Feynman gauge and use 
the on-mass-shell renormalization scheme. We compute the absorptive part of the Feynman diagram (which 
is free of the UV divergences) and then restore the real part using the appropriate dispersion relations. 
In the on-mass-shell scheme, the renormalized fermion self-energy vanishes like (j5q — m)^ which means 



that the second and third diagrams in Fig. 5.2c do not contribute to the cross section when the external quark 



legs are on-shell, pq m. The first graph in Fig. 5.2c describes the NLO corrections to the quark-photon 
vertex function: 

Am(9) = / (g') - ^^^z^-'?"' (5-62) 

where a^y = ^ {'J^ju — Ivl^i) while / iff') and g [ff') are the quark electromagnetic formfactors. At the 
lowest order A^°^ = 7^. 

The virtual lepton-quark cross section, ctYq, is obtained from the interference term between the virtual 
and the Born amplitude. The result can be written in terms of the electromagnetic formfactors as: 

^ = V'^'^^^ - { [1 + (1 - yf - 2A^^y1 / m + {Q') ] ■ (5.63) 



Taking into account the definition of the HM cross sections, o"*-^^ and a^^\ given by Eqs. (|p|), (pD and 
( p. 16 ), we find that corresponding virtual parts are: 



where /(^) (Q^) and 5(1) (Q^) 

are the NLO corrections to the electromagnetic formfactors. For the NLO 
HM cross sections, c|y and u^y , we use exactly the same notations as in Ref. ||38(| . 

In the on-mass-shell renormalization scheme, the renormalized vertex correction vanishes as the photon 
virtuality goes to zero, /^^^ (0) = 0. This is a consequence of the Ward identity and the fact that the real 
photon field (like the massive fermion one) is unrenormalized in first order QCD. To satisfy the condition 



/(I) (0) = automatically, we should use for f^^^ (g^) the dispersion relation with one subtraction. The 
second formfactor, g^^^ (g^) , has no singularities. For this reason, we use for g^^'^ (g^) the dispersion relation 
without subtractions: 
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TT 
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(5.65) 



Calculating the imaginary parts of the formafactors and restoring their real parts with the help of Eqs. ( |5.65| ) 
yields 
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(5.67) 



In Eqs. ( |5.66| ) and ( 5.67 ), Cp = {N^ — l)/(2iVc), where Nc is number of colors, and r is defined by Eq. ( ^.27| ). 
Taking into account that A^^)'^ = (a^ /aem) Cf ^q^eD' Fiqs. (|5.6(: , 5.67 ) reproduce the textbook 

QED results. 

It is now straightforward to obtain the virtual contribution to the longitudinal cross section. Combining 



Eqs. ( PI ) and (IOTD yields: 
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A Inr 
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(5.68) 



Comparing the above result with the corresponding one given by Eq. (39) in Ref. [^], we see that the HM 
expression for crly has opposite sign. Note also that this typo propagates into the final result for crj^^^ given 



byEq. (52) |3|]. 

Calculation of the bremsstrahlung contribution to the longitudinal cross section, c'"[^''(z, Q^), is also 
straightforward. We coincide with the HM result for (Tj^^(z,(5^) given by Eq. (49) in Ref. |3^. However 
there is one more misprint in the HM expression for a^"^ : the r.h.s of Eq. (52) should be multiplied by 
z. 

In the case of (7(^^(2;, Q^), the situation is slightly more complicated due to the need to take into account 
the IR singularities. One can see from Eq. (5.66) that /^^^ (Q'^) has an IR divergence which is regularized 
with the help of an infinitesimal gluon mass rUg. This singularity is cancelled when one adds the so-called soft 

contribution originating from the real gluon emission. For this purpose we introduce another infinitesimal 

(2) 

parameter {mg/m) ^ <C 1. The full bremsstrahlung contribution, <y\^-, can then be splitted into the 
soft and hard pieces as follows: 



aflMM')=e{z + 5z-l)a\%{z,Q% a'Zrd{^M')=0{l-z-5z)a\%{z,Q 



.(2) 



(2) 



oft 
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where 9(1 — z — 5z) is the Heaviside step function. The soft cross section should be calculated in the eikonal 
approximation, jTg ^ 0, taking into account the infinitesimal gluon mass nig. As a result, the sum of the 
virtual and soft contributions is IR finite: 



.(2) 
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a. 



' ^ 1 + 2A 

1 H , In r 

V1 + 4A 



+ 21nA-l-^^l±i5lnr 



(5.70) 



+ 



1 + 2A 



VI + 4A L 



Li2(r^) + 2Li2(-r) + -In^r - 21nr - Inr In A + 21nr ln(l + 4A) 
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Adding to the above expression the hard cross section cr^'^l^^^ defined by Eq. ( |5.69| ), we reproduce in the 
limit 5z ^ the full result for given by Eq. (51) in Ref. |38|. 



5.4 Appendix B: NLO Soft-Gluon Corrections to the Photon-Gluon Fu- 
sion 



This Appendix provides an overview of the NLO soft-gluon approximation for the photon-gluon fusion 
mechanism. We present the final results for the parton level cross sections to the next-to-leading logarithmic 
(NLL) accuracy. More details can be found in Refs. |36[| . 

To take into account the NLO contributions to the GF mechanism, one needs to calculate the virtual 
0{aemC(i) corrections to the Born process ( |5.39 ) and the real gluon emission: 



l*{q) + 9{kg) ^ Q{pq) + Q{Pq) + 9{Pg). 
The partonic invariants describing the single-particle inclusive (IPI) kinematics are 
s' = 2q ■ kg = s + Q"^ = CS', h = {kg - pqf - w? = (Ti, 



(5.71) 



S4 
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(5.72) 



where ( is defined by kg = (p and 54 measures the inelasticity of the reaction ( |5.7l| ). The corresponding 
IPI hadron level variables describing the reaction (5.1) 



are 



S' = 2q-p = S + Q^ 
84 = S' + T1 + U1, 



Ti = (p-pq) - m^ 
Ul = [q-pqf - m^. 



(5.73) 



The exact NLO calculations of the unpolarized heavy quark production in 717 |2^, |2^, j*g |2^, and gg 
p^ , |25| , |2^ , |27| collisions show that, near the partonic threshold, a strong logarithmic enhancement of the 
cross sections takes place in the collinear, pg^T — 0, and soft, pg 0, limits. This threshold (or soft-gluon) 
enhancement has universal nature in the perturbation theory and originates from incomplete cancellation 
of the soft and collinear singularities between the loop and the bremsstrahlung contributions. Large leading 
and next-to-leading threshold logarithms can be resummed to all orders of perturbative expansion using 
the appropriate evolution equations |28, 29, The analytic results for the resummed cross sections are 
ill-defined due to the Landau pole in the coupling strength a^. However, if one considers the obtained 
expressions as generating functionals of the perturbative theory and re-expands them at fixed order in Og, 
no divergences associated with the Landau pole are encountered. 



Soft-gluon resummation for the photon-gluon fusion has been performed in Ref. [31| and checked in 
Refs. [^, To NLL accuracy, the perturbative expansion for the partonic cross sections, d'^ak^g/dtidui 
(k = T, L, A, I), can be written in a factorized form as 

n=l 

with the Born level distributions B^'^'^ given by 



„/2 



^(5^^l,nl)=i?f°-(5',^l,Ul){5(/ + ^l+nl) + ^(^)"K(")(s^^l,nl)}, (5.74) 
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(5.78) 



Note that the functions K^"'\s' , ti, ui) in Eq. ( |5.74| ) originate from the colhnear and soft limits. Radiation 
of soft and colhnear gluons does not affect the transverse momentum of detected particles and therefore the 
azimuthal angle (f. For this reason, the functions K^'^\s' ,ti,ui) are the same for all helicity cross sections 
^k,g = T, L, A, I). At NLO, the soft-gluon corrections to NLL accuracy in the MS scheme are 
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where we use = fxp = fiR. In Eq. ( FiSOl) , Ca = and Cp = {N^ - l)/{2Nc), where is number of 

J+ivr} with (3 = \J\ — Am? / s. The single-particle inclusive 



colors, while Lp = (1 - 2mVs){ln[(l - /3)/(l + , 
"plus" distributions are defined by 



In' [si/m?) 


= lim < 


Si 





In' [si/rr? 

Si 



-9{si - e) + 



/ + 1 



In 



(5.80) 



For any sufficiently regular test function h{si), Eq. ( 5.80 ) gives 
In' [si/m"^) 



dsi h{si) 
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In Eq. ( [5.80| ), we have preserved the NLL terms for the scale-dependent logarithms too. Note also that the 
results ( |5.77D and ( 5.80| ) agree to NLL accuracy with the exact Oia^^c?^ calculations of the photon-gluon 
cross sections dT,g and dL,g given in Ref. [p^] . 

To investigate the scale dependence of the results ( 5.74 — SlSCj ), it is convenient to introduce for the 
fully inclusive (integrated over t\ and u\) cross sections, (T^^g {k = T, L, A, I), the dimensionless coefficient 



functions c 

.(1,0) 



(n,0 



(1,0) 



defined by Eq. ( pO]) . C oncerning the NLO scale-independent coefficient functions, only Cj, 



and c)^'^' are known exactly ||2^, 191]. As to the //-dependent coefficients, they can by calculated explicitly 
using the evolution equation: 
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din /i^ 
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(5.82) 



where z = Q'^/s', Cmin = z{l + 4A), ak,g{z,Q'^, fi) are the cross sections resummed to all orders in and 
PggiC) is the corresponding (resummed) Altarelli-Parisi gluon-gluon splitting function. Expanding Eq. ( |5.82| ) 
in tts, one can find ]31, [SSf] 
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where 62 = (llCyi — 2nj)/12 is the first coefficient of the /?(as)-function expansion and nj is the number of 
active quark flavors. The one-loop gluon splitting function is: 
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(5.84) 



With Eq. (5.83) in hand, it is possible to check the quality of the NLL approximation against exact 
answers. As shown in Ref. ]36], the soft-gluon corrections reproduce satisfactorily the threshold behavior 
of the available exact results for A 1. Since the gluon distribution function supports just the threshold 
region, the soft-gluon contribution dominates the photon-hadron cross sections (Tk^cF = T,L,A,I) at 
energies not so far from the production threshold and at relatively low virtuality < m^. 



5.5 Appendix C: Nonperturbative IC and Relevant Experimental Facts 



The most clean probe of the charm quark distribution function (both perturbative and nonperturbative) is 
the semi-inclusive deep inelastic lepton-proton scattering, Ip I'cX. To measure the nonperturbative IC 
contribution, one needs data on the charm production at sufficiently large Bjorken x. The only experiment 



which has investigated the large x domain is the European Muon Collaboration (EMC) [192] where the 



decay lepton spectra have been used to detect the produced charmed particles. In Ref. |193| , a re-analysis of 
the EMC data on (x, Q'^) have been performed using the NLO results for both GF and QS components. 
The analysis [193] shows that a nonperturbative intrinsic charm contribution to the proton wave function of 
the order of 1% is needed to fit the EMC data in the large x region. This value of the nonperturbative IC is 



consistent with the estimates based on the operator product expansion [169[. Note however that the EMC 
data are of limited statistics and, for this reason, more accurate measurements of charm leptoproduction at 
large x are necessary. 

It is also possible to extract useful information on the IC from diffractive dissociation processes such as 
p — > pJ/ip on a nuclear target. Comprehensive measurements of the pA — > J /ipX and ttA J/ipX cross 
sections have been performed in the fixed target experiments NA3 at CERN |194f[ and E886 at FNAL |195|]. 



According to the arguments presented in Refs. \19G, 197, 198], the IC contribution is predicted to be strongly 
shadowed in the above reactions that is in a complete agreement with the observed nuclear dependence of 
the high Feynman xp component of the J/ip hadroproduction. 

A non- vanishing five-quark Fock component \uudcc) leads to the production of open charm states such 
as Ac{cud) and D^{cd) with large Feynman xp. This may occur either through a coalescence of the valence 
and charm quarks which are moving with the same rapidity or via hadronization of the produced c and c. 



As shown in Refs. ]19£, 200], a model based on the nonperturbative intrinsic charm naturally explains the 
leading particle effect in the pp DX and pp — > KcX processes that has been observed at the ISR ] 201 ] 
and Fermilab |20|, |20^ ]. 

As to the high-Xi? hadroproduction of open bottom states like Ah{bud), corresponding cross sections are 
predicted to be suppressed as /ml ~ 1/10 in comparison with the case of charm production. Evidence 
for the forward A^ production in the pp collisions at the ISR energy was reported in Refs. ] ^04] , ^05| . 

Rare seven-quark fluctuations of the type \ uudcccc) in the proton wave function can lead to the production 
of two J/'ip [p06| or a double-charm baryon state at large xp and low pp- Double J/ip events with a high 



combined xp > 0.5 have been detected in the NA3 experiment |207| . An observation of the double-charmed 
baryon H+(3520) with mean (xp) ~ 0.33 has been reported by the SELEX collaboration at FNAL [p08[[ . 



Conclusion 



We conclude by summarizing our main observations and results. In this thesis, some spin effects in QCD 
and recurrence lattices with multi-site exchanges are investigated. Our conclusions concerning the critical 
phenomena in the recurrence spin models with multi-site exchanges can be formulated as follows. 

We have studied the Yang-Lee complex zeros of the partition function of the anti-ferromagnetic Potts 
model for biopolymer formulated on the recursive zigzag ladder with three-spin interaction. Using the 
dynamical system approach of multi-dimensional mapping, we have obtained the helix-coil pseudo-phase 
transitions in the thermodynamic limit. Taking into account a non-classical helix-stabilizing interaction, we 
describe the folding and unfolding processes for proteins and polypeptides. We have got also Arnold tongues 
with different winding numbers in the considered microscopic theory |7|. 



Using the dynamical system approach to solid and fluid '^He films with MSE on the recurrent lattices, 
we have obtained magnetization curves with plateaus (at m = 0, m = 1/2, m = 1/3 and m = 2/3) and 
one period doubling [|^, 82, |8^. We have taken into account two—, three—, four—, five—, and six— spin 
exchanges. Interaction parameters used in calculations were taken from experimental results. 

We have also considered (on the Bethe lattice) the model with cubic symmetry containing both linear 
and quadratic spin-spin interactions. An expression for the free energy per spin in the thermodynamic limit 
was found. Then we have identified the different thermodynamic phases of the system (disordered, partially 
ordered and completely ordered) in the ferromagnetic case {J > 0, K > 0). We have obtained the phase 
diagrams of the model which are found to be different for Q < 2 and Q > 2. The case Q < 2 contains 
three tricritical points. When Q > 2, one tricritical and one triple points there exist |81|. 

Next direction of our activity are spin effects in QCD. We consider the azimuthal dependence in charm 
leptoproduction as a probe of the IC content of the proton. Our analysis is based on the fact that the 
GF and QS components have strongly different cos 2(/?-distributions [^] . This fact follows from the NLO 
calculations of both parton level contributions. In the framework of the FFNS, we justify the most remarkable 
property of the hadron level azimuthal cos2(^ asymmetry: the combined GF+QS predictions for A^XjQ"^) 
are perturbatively and parametrically stable. The nonperturbative IC contribution (resulting from the five- 
quark \uudcc) component of the proton wave function) is practically invisible at low x, but affects essentially 
the GF predictions for the asymmetry at large Bjorken x. We conclude that measurements of the cos 2ip 
asymmetry at large x could directly probe the nonperturbative intrinsic charm p^]. 

Within the VFNS, charm density originates perturbatively from the 5 — > cc process and obeys the 
DGLAP evolution equation. Presently, charm densities are included practically in all the global sets of PDFs 
like CTEQ and MRST. Our analysis shows that these charm distribution functions reduce dramatically (by 
about 1/3) the GF predictions for A{x, Q"^) practically at all values of x. For this reason, the perturbative IC 
contribution can easily be measured using the azimuthal cos 2(^-distributions in charm leptoproduction . 



This dissertation is based on the material published in journals: Physics Letters A [79], Physical Review 
B []80| , Physica A |M| , International Journal of Modern Physics B |p2| , Proceedings of National Academy 



of Sciences of Armenia, Physics [^], Nuclear Physics B |8^, Physical Review D [^5|. Our main results 
are reported on conferences: Chaos and Supercomputers (Nor-Amberd, Armenia 2000), IX International 
Conference on Symmetry Methods in Physics (Yerevan, Armenia 2001), Stat. Mech. and Dynamic Systems 
(Nor-Amberd, Armenia 2003), "Fizika-2000" Young Scientist's Republic Conference (Yerevan, Armenia 
2004) , Selected Topics in Theoretical Physics (Tbilisi, Georgia 2005), and discussed on seminars in Yerevan 
Physics Institute and ICTP (Trieste, Italy). 
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